-----------------------------------------------------------------------------------------
      name:  <unnamed>
       log:  /Users/drdrtsai/Documents/Stata/Data/HopeNet/IPV_RCT/ipvrct_replication_2017
> 0408.log
  log type:  text
 opened on:   9 Apr 2017, 13:32:13

. 
. use ipvrct_replication_20170408.dta, replace

. 
. * summary statistics & covariate balance
. 
. sort visitdate 

. list visitdate in 1/5

     +-----------+
     | visitdate |
     |-----------|
  1. | 03jun2014 |
  2. | 03jun2014 |
  3. | 03jun2014 |
  4. | 03jun2014 |
  5. | 03jun2014 |
     +-----------+

. gsort -visitdate

. list visitdate in 1/5

     +-----------+
     | visitdate |
     |-----------|
  1. | 12aug2015 |
  2. | 22jul2015 |
  3. | 16jul2015 |
  4. | 07jul2015 |
  5. | 19may2015 |
     +-----------+

. 
. tab ipv_version, missing

        IPV |
 experiment |
  treatment |
        arm |      Freq.     Percent        Cum.
------------+-----------------------------------
1:DHS style |        454       34.03       34.03
    2:blame |        441       33.06       67.09
3:blameless |        439       32.91      100.00
------------+-----------------------------------
      Total |      1,334      100.00

. 
. tab female ipv_version, col chi2

+-------------------+
| Key               |
|-------------------|
|     frequency     |
| column percentage |
+-------------------+

     =1 if |
participan |
      t is |   IPV experiment treatment arm
    female | 1:DHS sty    2:blame  3:blamele |     Total
-----------+---------------------------------+----------
    0:male |       198        187        193 |       578 
           |     43.61      42.40      43.96 |     43.33 
-----------+---------------------------------+----------
  1:female |       256        254        246 |       756 
           |     56.39      57.60      56.04 |     56.67 
-----------+---------------------------------+----------
     Total |       454        441        439 |     1,334 
           |    100.00     100.00     100.00 |    100.00 

          Pearson chi2(2) =   0.2406   Pr = 0.887

. tabstat age, stat(n mean sd) by(ipv_version)

Summary for variables: age
     by categories of: ipv_version (IPV experiment treatment arm)

ipv_version |         N      mean        sd
------------+------------------------------
1:DHS style |       454  41.84361  17.10311
    2:blame |       441  41.47846   16.4636
3:blameless |       439  41.20046  17.22112
------------+------------------------------
      Total |      1334  41.51124  16.92314
-------------------------------------------

. ttest age if ipv_version==1 | ipv_version==2, by(ipv_version)

Two-sample t test with equal variances
------------------------------------------------------------------------------
   Group |     Obs        Mean    Std. Err.   Std. Dev.   [95% Conf. Interval]
---------+--------------------------------------------------------------------
1:DHS st |     454    41.84361    .8026887    17.10311    40.26616    43.42107
 2:blame |     441    41.47846    .7839809     16.4636    39.93765    43.01927
---------+--------------------------------------------------------------------
combined |     895    41.66369    .5609822    16.78265    40.56269    42.76468
---------+--------------------------------------------------------------------
    diff |            .3651543    1.122644               -1.838175    2.568483
------------------------------------------------------------------------------
    diff = mean(1:DHS st) - mean(2:blame)                         t =   0.3253
Ho: diff = 0                                     degrees of freedom =      893

    Ha: diff < 0                 Ha: diff != 0                 Ha: diff > 0
 Pr(T < t) = 0.6275         Pr(|T| > |t|) = 0.7451          Pr(T > t) = 0.3725

. ttest age if ipv_version==1 | ipv_version==3, by(ipv_version)

Two-sample t test with equal variances
------------------------------------------------------------------------------
   Group |     Obs        Mean    Std. Err.   Std. Dev.   [95% Conf. Interval]
---------+--------------------------------------------------------------------
1:DHS st |     454    41.84361    .8026887    17.10311    40.26616    43.42107
3:blamel |     439    41.20046    .8219193    17.22112    39.58506    42.81585
---------+--------------------------------------------------------------------
combined |     893    41.52744    .5740574    17.15462    40.40078     42.6541
---------+--------------------------------------------------------------------
    diff |            .6431568    1.148719               -1.611354    2.897667
------------------------------------------------------------------------------
    diff = mean(1:DHS st) - mean(3:blamel)                        t =   0.5599
Ho: diff = 0                                     degrees of freedom =      891

    Ha: diff < 0                 Ha: diff != 0                 Ha: diff > 0
 Pr(T < t) = 0.7122         Pr(|T| > |t|) = 0.5757          Pr(T > t) = 0.2878

. ttest age if ipv_version==2 | ipv_version==3, by(ipv_version)

Two-sample t test with equal variances
------------------------------------------------------------------------------
   Group |     Obs        Mean    Std. Err.   Std. Dev.   [95% Conf. Interval]
---------+--------------------------------------------------------------------
 2:blame |     441    41.47846    .7839809     16.4636    39.93765    43.01927
3:blamel |     439    41.20046    .8219193    17.22112    39.58506    42.81585
---------+--------------------------------------------------------------------
combined |     880    41.33977    .5675666    16.83675    40.22583    42.45372
---------+--------------------------------------------------------------------
    diff |            .2780025    1.135744               -1.951087    2.507092
------------------------------------------------------------------------------
    diff = mean(2:blame) - mean(3:blamel)                         t =   0.2448
Ho: diff = 0                                     degrees of freedom =      878

    Ha: diff < 0                 Ha: diff != 0                 Ha: diff > 0
 Pr(T < t) = 0.5967         Pr(|T| > |t|) = 0.8067          Pr(T > t) = 0.4033

. tab edcat ipv_version, col chi2

+-------------------+
| Key               |
|-------------------|
|     frequency     |
| column percentage |
+-------------------+

        Participant's |
          educational |   IPV experiment treatment arm
           attainment | 1:DHS sty    2:blame  3:blamele |     Total
----------------------+---------------------------------+----------
               0:none |        83         89         73 |       245 
                      |     18.28      20.18      16.63 |     18.37 
----------------------+---------------------------------+----------
1:some primary (P1-P6 |       140        136        148 |       424 
                      |     30.84      30.84      33.71 |     31.78 
----------------------+---------------------------------+----------
2:completed primary ( |       110        108         98 |       316 
                      |     24.23      24.49      22.32 |     23.69 
----------------------+---------------------------------+----------
3:primary+ (S1-S6,voc |       121        108        120 |       349 
                      |     26.65      24.49      27.33 |     26.16 
----------------------+---------------------------------+----------
                Total |       454        441        439 |     1,334 
                      |    100.00     100.00     100.00 |    100.00 

          Pearson chi2(6) =   3.5452   Pr = 0.738

. tab mstatus ipv_version, col chi2

+-------------------+
| Key               |
|-------------------|
|     frequency     |
| column percentage |
+-------------------+

Participant's marital |   IPV experiment treatment arm
               status | 1:DHS sty    2:blame  3:blamele |     Total
----------------------+---------------------------------+----------
1:married/cohabitatin |       299        285        301 |       885 
                      |     66.00      64.63      68.56 |     66.39 
----------------------+---------------------------------+----------
2:separated,divorced, |        95         95         80 |       270 
                      |     20.97      21.54      18.22 |     20.26 
----------------------+---------------------------------+----------
3:single,never marrie |        59         61         58 |       178 
                      |     13.02      13.83      13.21 |     13.35 
----------------------+---------------------------------+----------
                Total |       453        441        439 |     1,333 
                      |    100.00     100.00     100.00 |    100.00 

          Pearson chi2(4) =   2.0187   Pr = 0.732

. tab vid ipv_version, col chi2

+-------------------+
| Key               |
|-------------------|
|     frequency     |
| column percentage |
+-------------------+

           |   IPV experiment treatment arm
village ID | 1:DHS sty    2:blame  3:blamele |     Total
-----------+---------------------------------+----------
         1 |        61         67         63 |       191 
           |     13.44      15.19      14.35 |     14.32 
-----------+---------------------------------+----------
         2 |        64         64         61 |       189 
           |     14.10      14.51      13.90 |     14.17 
-----------+---------------------------------+----------
         3 |        62         60         58 |       180 
           |     13.66      13.61      13.21 |     13.49 
-----------+---------------------------------+----------
         4 |        61         55         57 |       173 
           |     13.44      12.47      12.98 |     12.97 
-----------+---------------------------------+----------
         5 |        40         36         39 |       115 
           |      8.81       8.16       8.88 |      8.62 
-----------+---------------------------------+----------
         6 |        60         61         65 |       186 
           |     13.22      13.83      14.81 |     13.94 
-----------+---------------------------------+----------
         7 |        41         37         38 |       116 
           |      9.03       8.39       8.66 |      8.70 
-----------+---------------------------------+----------
         8 |        65         61         58 |       184 
           |     14.32      13.83      13.21 |     13.79 
-----------+---------------------------------+----------
     Total |       454        441        439 |     1,334 
           |    100.00     100.00     100.00 |    100.00 

         Pearson chi2(14) =   1.6209   Pr = 1.000

. mlogit ipv_version female age i.mstatus i.vid

Iteration 0:   log likelihood = -1464.3215  
Iteration 1:   log likelihood = -1462.3157  
Iteration 2:   log likelihood =  -1462.315  
Iteration 3:   log likelihood =  -1462.315  

Multinomial logistic regression                 Number of obs     =      1,333
                                                LR chi2(22)       =       4.01
                                                Prob > chi2       =     1.0000
Log likelihood =  -1462.315                     Pseudo R2         =     0.0014

----------------------------------------------------------------------------------------
           ipv_version |      Coef.   Std. Err.      z    P>|z|     [95% Conf. Interval]
-----------------------+----------------------------------------------------------------
1_DHS_style            |  (base outcome)
-----------------------+----------------------------------------------------------------
2_blame                |
                female |   .0496486    .142202     0.35   0.727    -.2290622    .3283593
                   age |  -.0013125   .0045093    -0.29   0.771    -.0101506    .0075256
                       |
               mstatus |
2:separated,divorce~w  |   .0463873   .1841883     0.25   0.801    -.3146151    .4073897
3:single,never marr..  |   .0652171   .2142452     0.30   0.761    -.3546958    .4851299
                       |
                   vid |
                    2  |   -.096416   .2504302    -0.39   0.700    -.5872502    .3944183
                    3  |  -.1299654   .2535122    -0.51   0.608    -.6268401    .3669093
                    4  |   -.197665   .2569376    -0.77   0.442    -.7012534    .3059235
                    5  |  -.2027726   .2901902    -0.70   0.485    -.7715349    .3659898
                    6  |  -.0868748   .2554074    -0.34   0.734    -.5874641    .4137145
                    7  |  -.1692207   .2888228    -0.59   0.558     -.735303    .3968615
                    8  |   -.157082   .2520574    -0.62   0.533    -.6511056    .3369415
                       |
                 _cons |   .1038319   .2749203     0.38   0.706    -.4350019    .6426657
-----------------------+----------------------------------------------------------------
3_blameless            |
                female |   .0148433   .1417944     0.10   0.917    -.2630686    .2927552
                   age |  -.0015136   .0045481    -0.33   0.739    -.0104277    .0074004
                       |
               mstatus |
2:separated,divorce~w  |  -.1733989   .1891145    -0.92   0.359    -.5440564    .1972587
3:single,never marr..  |  -.0625104   .2159436    -0.29   0.772    -.4857521    .3607313
                       |
                   vid |
                    2  |  -.0705367   .2539058    -0.28   0.781     -.568183    .4271097
                    3  |  -.1028967   .2565987    -0.40   0.688     -.605821    .4000276
                    4  |  -.0958019   .2576092    -0.37   0.710    -.6007067    .4091028
                    5  |  -.0586112   .2882275    -0.20   0.839    -.6235267    .5063043
                    6  |   .0558757    .255321     0.22   0.827    -.4445443    .5562956
                    7  |  -.0790716   .2893286    -0.27   0.785    -.6461452    .4880021
                    8  |  -.1665469   .2556431    -0.65   0.515    -.6675981    .3345044
                       |
                 _cons |   .1288828   .2776412     0.46   0.642     -.415284    .6730497
----------------------------------------------------------------------------------------

. test female age 2.mstatus 3.mstatus /*
> */ 2.vid 3.vid 4.vid 5.vid 6.vid 7.vid 8.vid

 ( 1)  [1_DHS_style]o.female = 0
 ( 2)  [2_blame]female = 0
 ( 3)  [3_blameless]female = 0
 ( 4)  [1_DHS_style]o.age = 0
 ( 5)  [2_blame]age = 0
 ( 6)  [3_blameless]age = 0
 ( 7)  [1_DHS_style]2o.mstatus = 0
 ( 8)  [2_blame]2.mstatus = 0
 ( 9)  [3_blameless]2.mstatus = 0
 (10)  [1_DHS_style]3o.mstatus = 0
 (11)  [2_blame]3.mstatus = 0
 (12)  [3_blameless]3.mstatus = 0
 (13)  [1_DHS_style]2o.vid = 0
 (14)  [2_blame]2.vid = 0
 (15)  [3_blameless]2.vid = 0
 (16)  [1_DHS_style]3o.vid = 0
 (17)  [2_blame]3.vid = 0
 (18)  [3_blameless]3.vid = 0
 (19)  [1_DHS_style]4o.vid = 0
 (20)  [2_blame]4.vid = 0
 (21)  [3_blameless]4.vid = 0
 (22)  [1_DHS_style]5o.vid = 0
 (23)  [2_blame]5.vid = 0
 (24)  [3_blameless]5.vid = 0
 (25)  [1_DHS_style]6o.vid = 0
 (26)  [2_blame]6.vid = 0
 (27)  [3_blameless]6.vid = 0
 (28)  [1_DHS_style]7o.vid = 0
 (29)  [2_blame]7.vid = 0
 (30)  [3_blameless]7.vid = 0
 (31)  [1_DHS_style]8o.vid = 0
 (32)  [2_blame]8.vid = 0
 (33)  [3_blameless]8.vid = 0
       Constraint 1 dropped
       Constraint 4 dropped
       Constraint 7 dropped
       Constraint 10 dropped
       Constraint 13 dropped
       Constraint 16 dropped
       Constraint 19 dropped
       Constraint 22 dropped
       Constraint 25 dropped
       Constraint 28 dropped
       Constraint 31 dropped

           chi2( 22) =    3.98
         Prob > chi2 =    1.0000

. 
. * correlate norms & beliefs
. 
. spearman belief_tot norm_tot

 Number of obs =    1334
Spearman's rho =       0.4988

Test of Ho: belief_tot and norm_tot are independent
    Prob > |t| =       0.0000

. 
. * first examine frequencies across arms
. 
. foreach x of numlist 1/5 {
  2. tab belief`x' ipv_version, col
  3. tab norm`x' ipv_version, col
  4. }

+-------------------+
| Key               |
|-------------------|
|     frequency     |
| column percentage |
+-------------------+

     Abuse |
 justified |
   if wife |
  goes out |
   without |   IPV experiment treatment arm
permission | 1:DHS sty    2:blame  3:blamele |     Total
-----------+---------------------------------+----------
         0 |       351        260        347 |       958 
           |     77.31      58.96      79.04 |     71.81 
-----------+---------------------------------+----------
         1 |       103        181         92 |       376 
           |     22.69      41.04      20.96 |     28.19 
-----------+---------------------------------+----------
     Total |       454        441        439 |     1,334 
           |    100.00     100.00     100.00 |    100.00 


+-------------------+
| Key               |
|-------------------|
|     frequency     |
| column percentage |
+-------------------+

>1/2 would |
     agree |
     abuse |
 justified |
   if wife |
  goes out |
   without |   IPV experiment treatment arm
permission | 1:DHS sty    2:blame  3:blamele |     Total
-----------+---------------------------------+----------
         0 |       299        235        317 |       851 
           |     65.86      53.29      72.21 |     63.79 
-----------+---------------------------------+----------
         1 |       155        206        122 |       483 
           |     34.14      46.71      27.79 |     36.21 
-----------+---------------------------------+----------
     Total |       454        441        439 |     1,334 
           |    100.00     100.00     100.00 |    100.00 


+-------------------+
| Key               |
|-------------------|
|     frequency     |
| column percentage |
+-------------------+

     Abuse |
 justified |
   if wife |
  neglects |   IPV experiment treatment arm
  children | 1:DHS sty    2:blame  3:blamele |     Total
-----------+---------------------------------+----------
         0 |       358        170        384 |       912 
           |     78.85      38.55      87.47 |     68.37 
-----------+---------------------------------+----------
         1 |        96        271         55 |       422 
           |     21.15      61.45      12.53 |     31.63 
-----------+---------------------------------+----------
     Total |       454        441        439 |     1,334 
           |    100.00     100.00     100.00 |    100.00 


+-------------------+
| Key               |
|-------------------|
|     frequency     |
| column percentage |
+-------------------+

>1/2 would |
     agree |
     abuse |
 justified |
   if wife |
  neglects |   IPV experiment treatment arm
  children | 1:DHS sty    2:blame  3:blamele |     Total
-----------+---------------------------------+----------
         0 |       336        200        364 |       900 
           |     74.01      45.35      82.92 |     67.47 
-----------+---------------------------------+----------
         1 |       118        241         75 |       434 
           |     25.99      54.65      17.08 |     32.53 
-----------+---------------------------------+----------
     Total |       454        441        439 |     1,334 
           |    100.00     100.00     100.00 |    100.00 


+-------------------+
| Key               |
|-------------------|
|     frequency     |
| column percentage |
+-------------------+

     Abuse |
 justified |
   if wife |
 argues in |   IPV experiment treatment arm
    public | 1:DHS sty    2:blame  3:blamele |     Total
-----------+---------------------------------+----------
         0 |       304        147        355 |       806 
           |     66.96      33.33      80.87 |     60.42 
-----------+---------------------------------+----------
         1 |       150        294         84 |       528 
           |     33.04      66.67      19.13 |     39.58 
-----------+---------------------------------+----------
     Total |       454        441        439 |     1,334 
           |    100.00     100.00     100.00 |    100.00 


+-------------------+
| Key               |
|-------------------|
|     frequency     |
| column percentage |
+-------------------+

>1/2 would |
     agree |
     abuse |
 justified |
   if wife |
 argues in |   IPV experiment treatment arm
    public | 1:DHS sty    2:blame  3:blamele |     Total
-----------+---------------------------------+----------
         0 |       271        179        307 |       757 
           |     59.69      40.59      69.93 |     56.75 
-----------+---------------------------------+----------
         1 |       183        262        132 |       577 
           |     40.31      59.41      30.07 |     43.25 
-----------+---------------------------------+----------
     Total |       454        441        439 |     1,334 
           |    100.00     100.00     100.00 |    100.00 


+-------------------+
| Key               |
|-------------------|
|     frequency     |
| column percentage |
+-------------------+

     Abuse |
 justified |
   if wife |
   refuses |
intercours |   IPV experiment treatment arm
         e | 1:DHS sty    2:blame  3:blamele |     Total
-----------+---------------------------------+----------
         0 |       378        275        376 |     1,029 
           |     83.26      62.36      85.65 |     77.14 
-----------+---------------------------------+----------
         1 |        76        166         63 |       305 
           |     16.74      37.64      14.35 |     22.86 
-----------+---------------------------------+----------
     Total |       454        441        439 |     1,334 
           |    100.00     100.00     100.00 |    100.00 


+-------------------+
| Key               |
|-------------------|
|     frequency     |
| column percentage |
+-------------------+

>1/2 would |
     agree |
     abuse |
 justified |
   if wife |
   refuses |
intercours |   IPV experiment treatment arm
         e | 1:DHS sty    2:blame  3:blamele |     Total
-----------+---------------------------------+----------
         0 |       311        260        315 |       886 
           |     68.50      58.96      71.75 |     66.42 
-----------+---------------------------------+----------
         1 |       143        181        124 |       448 
           |     31.50      41.04      28.25 |     33.58 
-----------+---------------------------------+----------
     Total |       454        441        439 |     1,334 
           |    100.00     100.00     100.00 |    100.00 


+-------------------+
| Key               |
|-------------------|
|     frequency     |
| column percentage |
+-------------------+

     Abuse |
 justified |
   if wife |
  does not |
   prepare |
   food on |   IPV experiment treatment arm
      time | 1:DHS sty    2:blame  3:blamele |     Total
-----------+---------------------------------+----------
         0 |       307        145        393 |       845 
           |     67.62      32.88      89.52 |     63.34 
-----------+---------------------------------+----------
         1 |       147        296         46 |       489 
           |     32.38      67.12      10.48 |     36.66 
-----------+---------------------------------+----------
     Total |       454        441        439 |     1,334 
           |    100.00     100.00     100.00 |    100.00 


+-------------------+
| Key               |
|-------------------|
|     frequency     |
| column percentage |
+-------------------+

>1/2 would |
     agree |
     abuse |
 justified |
   if wife |
  does not |
   prepare |
   food on |   IPV experiment treatment arm
      time | 1:DHS sty    2:blame  3:blamele |     Total
-----------+---------------------------------+----------
         0 |       320        187        371 |       878 
           |     70.48      42.40      84.51 |     65.82 
-----------+---------------------------------+----------
         1 |       134        254         68 |       456 
           |     29.52      57.60      15.49 |     34.18 
-----------+---------------------------------+----------
     Total |       454        441        439 |     1,334 
           |    100.00     100.00     100.00 |    100.00 


. 
. * examine means and totals (assuming DNK/refuse=0)
. 
. tabstat belief_tot, stat(n mean sd variance) by(ipv_version)

Summary for variables: belief_tot
     by categories of: ipv_version (IPV experiment treatment arm)

ipv_version |         N      mean        sd  variance
------------+----------------------------------------
1:DHS style |       454  1.259912  1.581718  2.501833
    2:blame |       441  2.739229  1.810959  3.279571
3:blameless |       439  .7744875  1.188486  1.412498
------------+----------------------------------------
      Total |      1334  1.589205  1.757821  3.089936
-----------------------------------------------------

. tabstat norm_tot, stat(n mean sd variance) by(ipv_version)

Summary for variables: norm_tot
     by categories of: ipv_version (IPV experiment treatment arm)

ipv_version |         N      mean        sd  variance
------------+----------------------------------------
1:DHS style |       454  1.614537  1.680655  2.824601
    2:blame |       441  2.594104  1.917826  3.678056
3:blameless |       439  1.186788  1.487147  2.211606
------------+----------------------------------------
      Total |      1334  1.797601   1.80134  3.244825
-----------------------------------------------------

. 
. ***** primary analysis: beliefs
. 
. * establish violation of proportional odds assumption
. 
. ologit belief_tot i.ipv_version i.female i.vid, vce(cluster vid) or

Iteration 0:   log pseudolikelihood = -2121.6536  
Iteration 1:   log pseudolikelihood = -1960.4571  
Iteration 2:   log pseudolikelihood = -1957.4203  
Iteration 3:   log pseudolikelihood = -1957.4123  
Iteration 4:   log pseudolikelihood = -1957.4123  

Ordered logistic regression                     Number of obs     =      1,334
                                                Wald chi2(2)      =          .
                                                Prob > chi2       =          .
Log pseudolikelihood = -1957.4123               Pseudo R2         =     0.0774

                                    (Std. Err. adjusted for 8 clusters in vid)
------------------------------------------------------------------------------
             |               Robust
  belief_tot | Odds Ratio   Std. Err.      z    P>|z|     [95% Conf. Interval]
-------------+----------------------------------------------------------------
 ipv_version |
    2:blame  |   4.997398   .3862802    20.81   0.000     4.294862    5.814851
3:blameless  |   .5940158   .1069739    -2.89   0.004     .4173578    .8454491
             |
      female |
   1:female  |   1.818801   .1247642     8.72   0.000     1.589994    2.080535
             |
         vid |
          2  |   .8899391   .0070165   -14.79   0.000     .8762927    .9037979
          3  |   1.218525   .0085099    28.30   0.000     1.201959    1.235319
          4  |   1.164641   .0076671    23.15   0.000      1.14971    1.179766
          5  |   .8901065   .0098232   -10.55   0.000     .8710602    .9095694
          6  |   .8006715   .0090396   -19.69   0.000     .7831488    .8185862
          7  |   1.109825   .0041878    27.61   0.000     1.101647    1.118063
          8  |   1.153039    .004087    40.17   0.000     1.145056    1.161077
-------------+----------------------------------------------------------------
       /cut1 |   .3136421   .0565795                      .2027483    .4245359
       /cut2 |   1.026097   .0747427                      .8796038     1.17259
       /cut3 |   1.660757   .0670482                      1.529345    1.792169
       /cut4 |   2.350957   .1066441                      2.141938    2.559975
       /cut5 |   3.370899   .0864887                      3.201384    3.540413
------------------------------------------------------------------------------

. 
. tab ipv_version, generate(ipv_version)

        IPV |
 experiment |
  treatment |
        arm |      Freq.     Percent        Cum.
------------+-----------------------------------
1:DHS style |        454       34.03       34.03
    2:blame |        441       33.06       67.09
3:blameless |        439       32.91      100.00
------------+-----------------------------------
      Total |      1,334      100.00

. tab vid, generate(vid)

 village ID |      Freq.     Percent        Cum.
------------+-----------------------------------
          1 |        191       14.32       14.32
          2 |        189       14.17       28.49
          3 |        180       13.49       41.98
          4 |        173       12.97       54.95
          5 |        115        8.62       63.57
          6 |        186       13.94       77.51
          7 |        116        8.70       86.21
          8 |        184       13.79      100.00
------------+-----------------------------------
      Total |      1,334      100.00

. quietly ologit belief_tot ipv_version2 ipv_version3 female vid2 vid3 vid4 vid5 vid6 vid
> 7 vid8, vce(cluster vid)

. brant

Brant Test of Parallel Regression Assumption

    Variable |      chi2   p>chi2    df
-------------+--------------------------
         All |     73.05    0.001    40
-------------+--------------------------
ipv_version2 |     11.32    0.023     4
ipv_version3 |     20.79    0.000     4
      female |      2.04    0.729     4
        vid2 |      3.30    0.509     4
        vid3 |      3.16    0.531     4
        vid4 |     10.06    0.039     4
        vid5 |      4.42    0.352     4
        vid6 |      2.11    0.715     4
        vid7 |      1.25    0.869     4
        vid8 |      4.90    0.298     4
----------------------------------------

A significant test statistic provides evidence that the parallel
regression assumption has been violated.

. 
. * switch to partial proportional odds
. 
. gologit2 belief_tot ipv_version2 ipv_version3 female vid2 vid3 vid4 vid5 vid6 vid7 vid8
> , /*
> */ robust cluster(vid) pl(female vid2 vid3 vid4 vid5 vid6 vid7 vid8) or

Generalized Ordered Logit Estimates             Number of obs     =      1,334
                                                Wald chi2(2)      =      61.46
                                                Prob > chi2       =     0.0000
Log pseudolikelihood = -1936.1223               Pseudo R2         =     0.0874

 ( 1)  [0]female - [1]female = 0
 ( 2)  [0]vid2 - [1]vid2 = 0
 ( 3)  [0]vid3 - [1]vid3 = 0
 ( 4)  [0]vid4 - [1]vid4 = 0
 ( 5)  [0]vid5 - [1]vid5 = 0
 ( 6)  [0]vid6 - [1]vid6 = 0
 ( 7)  [0]vid7 - [1]vid7 = 0
 ( 8)  [0]vid8 - [1]vid8 = 0
 ( 9)  [1]female - [2]female = 0
 (10)  [1]vid2 - [2]vid2 = 0
 (11)  [1]vid3 - [2]vid3 = 0
 (12)  [1]vid4 - [2]vid4 = 0
 (13)  [1]vid5 - [2]vid5 = 0
 (14)  [1]vid6 - [2]vid6 = 0
 (15)  [1]vid7 - [2]vid7 = 0
 (16)  [1]vid8 - [2]vid8 = 0
 (17)  [2]female - [3]female = 0
 (18)  [2]vid2 - [3]vid2 = 0
 (19)  [2]vid3 - [3]vid3 = 0
 (20)  [2]vid4 - [3]vid4 = 0
 (21)  [2]vid5 - [3]vid5 = 0
 (22)  [2]vid6 - [3]vid6 = 0
 (23)  [2]vid7 - [3]vid7 = 0
 (24)  [2]vid8 - [3]vid8 = 0
 (25)  [3]female - [4]female = 0
 (26)  [3]vid2 - [4]vid2 = 0
 (27)  [3]vid3 - [4]vid3 = 0
 (28)  [3]vid4 - [4]vid4 = 0
 (29)  [3]vid5 - [4]vid5 = 0
 (30)  [3]vid6 - [4]vid6 = 0
 (31)  [3]vid7 - [4]vid7 = 0
 (32)  [3]vid8 - [4]vid8 = 0
                                    (Std. Err. adjusted for 8 clusters in vid)
------------------------------------------------------------------------------
             |               Robust
  belief_tot | Odds Ratio   Std. Err.      z    P>|z|     [95% Conf. Interval]
-------------+----------------------------------------------------------------
0            |
ipv_version2 |   3.870481   .4706286    11.13   0.000     3.049746    4.912089
ipv_version3 |   .6955832   .1414289    -1.79   0.074     .4669595    1.036141
      female |    1.78676   .1229564     8.43   0.000     1.561315    2.044758
        vid2 |   .8942266   .0079062   -12.64   0.000     .8788641    .9098576
        vid3 |   1.225832   .0092934    26.86   0.000     1.207752    1.244182
        vid4 |   1.168219   .0072855    24.93   0.000     1.154027    1.182586
        vid5 |   .8930522   .0092755   -10.89   0.000     .8750564    .9114182
        vid6 |   .7848355   .0094247   -20.18   0.000     .7665791    .8035267
        vid7 |   1.085399   .0045401    19.59   0.000     1.076537    1.094334
        vid8 |    1.12887   .0058673    23.32   0.000     1.117429    1.140429
       _cons |   .7144905    .058626    -4.10   0.000     .6083492    .8391507
-------------+----------------------------------------------------------------
1            |
ipv_version2 |   4.562102   .3983982    17.38   0.000     3.844425    5.413754
ipv_version3 |   .4129059   .0810134    -4.51   0.000     .2810883      .60654
      female |    1.78676   .1229564     8.43   0.000     1.561315    2.044758
        vid2 |   .8942266   .0079062   -12.64   0.000     .8788641    .9098576
        vid3 |   1.225832   .0092934    26.86   0.000     1.207752    1.244182
        vid4 |   1.168219   .0072855    24.93   0.000     1.154027    1.182586
        vid5 |   .8930522   .0092755   -10.89   0.000     .8750564    .9114182
        vid6 |   .7848355   .0094247   -20.18   0.000     .7665791    .8035267
        vid7 |   1.085399   .0045401    19.59   0.000     1.076537    1.094334
        vid8 |    1.12887   .0058673    23.32   0.000     1.117429    1.140429
       _cons |   .3945577   .0218598   -16.79   0.000     .3539575    .4398148
-------------+----------------------------------------------------------------
2            |
ipv_version2 |   5.738971   .7740374    12.95   0.000     4.405845    7.475477
ipv_version3 |   .4497126   .1107531    -3.24   0.001     .2775269    .7287271
      female |    1.78676   .1229564     8.43   0.000     1.561315    2.044758
        vid2 |   .8942266   .0079062   -12.64   0.000     .8788641    .9098576
        vid3 |   1.225832   .0092934    26.86   0.000     1.207752    1.244182
        vid4 |   1.168219   .0072855    24.93   0.000     1.154027    1.182586
        vid5 |   .8930522   .0092755   -10.89   0.000     .8750564    .9114182
        vid6 |   .7848355   .0094247   -20.18   0.000     .7665791    .8035267
        vid7 |   1.085399   .0045401    19.59   0.000     1.076537    1.094334
        vid8 |    1.12887   .0058673    23.32   0.000     1.117429    1.140429
       _cons |   .1856349   .0190222   -16.43   0.000     .1518575    .2269255
-------------+----------------------------------------------------------------
3            |
ipv_version2 |   5.184952   .9646944     8.85   0.000     3.600586    7.466489
ipv_version3 |   .3642585   .1154273    -3.19   0.001     .1957395    .6778615
      female |    1.78676   .1229564     8.43   0.000     1.561315    2.044758
        vid2 |   .8942266   .0079062   -12.64   0.000     .8788641    .9098576
        vid3 |   1.225832   .0092934    26.86   0.000     1.207752    1.244182
        vid4 |   1.168219   .0072855    24.93   0.000     1.154027    1.182586
        vid5 |   .8930522   .0092755   -10.89   0.000     .8750564    .9114182
        vid6 |   .7848355   .0094247   -20.18   0.000     .7665791    .8035267
        vid7 |   1.085399   .0045401    19.59   0.000     1.076537    1.094334
        vid8 |    1.12887   .0058673    23.32   0.000     1.117429    1.140429
       _cons |   .1004941   .0167174   -13.81   0.000      .072534    .1392322
-------------+----------------------------------------------------------------
4            |
ipv_version2 |   3.908104   .9965898     5.35   0.000     2.370853    6.442101
ipv_version3 |   .2499627   .1027657    -3.37   0.001     .1116672    .5595315
      female |    1.78676   .1229564     8.43   0.000     1.561315    2.044758
        vid2 |   .8942266   .0079062   -12.64   0.000     .8788641    .9098576
        vid3 |   1.225832   .0092934    26.86   0.000     1.207752    1.244182
        vid4 |   1.168219   .0072855    24.93   0.000     1.154027    1.182586
        vid5 |   .8930522   .0092755   -10.89   0.000     .8750564    .9114182
        vid6 |   .7848355   .0094247   -20.18   0.000     .7665791    .8035267
        vid7 |   1.085399   .0045401    19.59   0.000     1.076537    1.094334
        vid8 |    1.12887   .0058673    23.32   0.000     1.117429    1.140429
       _cons |   .0447634   .0079204   -17.56   0.000     .0316456     .063319
------------------------------------------------------------------------------

. 
. ***** primary analysis: norms
. 
. * establish violation of proportional odds assumption
. 
. ologit norm_tot i.ipv_version i.female i.vid, vce(cluster vid) or

Iteration 0:   log pseudolikelihood = -2216.2268  
Iteration 1:   log pseudolikelihood = -2132.2308  
Iteration 2:   log pseudolikelihood = -2131.6439  
Iteration 3:   log pseudolikelihood = -2131.6436  
Iteration 4:   log pseudolikelihood = -2131.6436  

Ordered logistic regression                     Number of obs     =      1,334
                                                Wald chi2(2)      =          .
                                                Prob > chi2       =          .
Log pseudolikelihood = -2131.6436               Pseudo R2         =     0.0382

                                    (Std. Err. adjusted for 8 clusters in vid)
------------------------------------------------------------------------------
             |               Robust
    norm_tot | Odds Ratio   Std. Err.      z    P>|z|     [95% Conf. Interval]
-------------+----------------------------------------------------------------
 ipv_version |
    2:blame  |   2.624496   .2470293    10.25   0.000     2.182363    3.156202
3:blameless  |   .6255025   .0505514    -5.81   0.000     .5338721    .7328597
             |
      female |
   1:female  |   1.688083   .1987179     4.45   0.000     1.340269    2.126159
             |
         vid |
          2  |   .8432295   .0082142   -17.50   0.000     .8272826    .8594838
          3  |   1.145008   .0027029    57.36   0.000     1.139723    1.150318
          4  |   1.013473   .0044657     3.04   0.002     1.004759    1.022264
          5  |   .8638101   .0092601   -13.66   0.000       .84585    .8821515
          6  |   1.089616   .0109063     8.57   0.000     1.068448    1.111203
          7  |   .8618546   .0056548   -22.66   0.000     .8508423    .8730093
          8  |   1.249965   .0099058    28.15   0.000       1.2307    1.269532
-------------+----------------------------------------------------------------
       /cut1 |  -.0939561   .0725578                     -.2361668    .0482546
       /cut2 |   .5329873   .0474713                      .4399452    .6260293
       /cut3 |   1.124322    .061553                       1.00368    1.244963
       /cut4 |   1.768703   .0757379                      1.620259    1.917146
       /cut5 |   2.656126   .1102018                      2.440134    2.872117
------------------------------------------------------------------------------

. 
. quietly ologit norm_tot ipv_version2 ipv_version3 female vid2 vid3 vid4 vid5 vid6 vid7 
> vid8, vce(cluster vid)

. brant

Brant Test of Parallel Regression Assumption

    Variable |      chi2   p>chi2    df
-------------+--------------------------
         All |     43.64    0.319    40
-------------+--------------------------
ipv_version2 |     10.74    0.030     4
ipv_version3 |      2.16    0.706     4
      female |      1.91    0.753     4
        vid2 |      2.33    0.676     4
        vid3 |      2.28    0.685     4
        vid4 |      0.19    0.996     4
        vid5 |      1.73    0.785     4
        vid6 |      2.00    0.735     4
        vid7 |      0.52    0.971     4
        vid8 |      8.81    0.066     4
----------------------------------------

A significant test statistic provides evidence that the parallel
regression assumption has been violated.

. 
. * switch to partial proportional odds
. 
. gologit2 norm_tot ipv_version2 ipv_version3 female vid2 vid3 vid4 vid5 vid6 vid7 vid8, 
> /*
> */ robust cluster(vid) pl(female vid2 vid3 vid4 vid5 vid6 vid7 vid8) or

Generalized Ordered Logit Estimates             Number of obs     =      1,334
                                                Wald chi2(2)      =      45.61
                                                Prob > chi2       =     0.0000
Log pseudolikelihood = -2121.8375               Pseudo R2         =     0.0426

 ( 1)  [0]female - [1]female = 0
 ( 2)  [0]vid2 - [1]vid2 = 0
 ( 3)  [0]vid3 - [1]vid3 = 0
 ( 4)  [0]vid4 - [1]vid4 = 0
 ( 5)  [0]vid5 - [1]vid5 = 0
 ( 6)  [0]vid6 - [1]vid6 = 0
 ( 7)  [0]vid7 - [1]vid7 = 0
 ( 8)  [0]vid8 - [1]vid8 = 0
 ( 9)  [1]female - [2]female = 0
 (10)  [1]vid2 - [2]vid2 = 0
 (11)  [1]vid3 - [2]vid3 = 0
 (12)  [1]vid4 - [2]vid4 = 0
 (13)  [1]vid5 - [2]vid5 = 0
 (14)  [1]vid6 - [2]vid6 = 0
 (15)  [1]vid7 - [2]vid7 = 0
 (16)  [1]vid8 - [2]vid8 = 0
 (17)  [2]female - [3]female = 0
 (18)  [2]vid2 - [3]vid2 = 0
 (19)  [2]vid3 - [3]vid3 = 0
 (20)  [2]vid4 - [3]vid4 = 0
 (21)  [2]vid5 - [3]vid5 = 0
 (22)  [2]vid6 - [3]vid6 = 0
 (23)  [2]vid7 - [3]vid7 = 0
 (24)  [2]vid8 - [3]vid8 = 0
 (25)  [3]female - [4]female = 0
 (26)  [3]vid2 - [4]vid2 = 0
 (27)  [3]vid3 - [4]vid3 = 0
 (28)  [3]vid4 - [4]vid4 = 0
 (29)  [3]vid5 - [4]vid5 = 0
 (30)  [3]vid6 - [4]vid6 = 0
 (31)  [3]vid7 - [4]vid7 = 0
 (32)  [3]vid8 - [4]vid8 = 0
                                    (Std. Err. adjusted for 8 clusters in vid)
------------------------------------------------------------------------------
             |               Robust
    norm_tot | Odds Ratio   Std. Err.      z    P>|z|     [95% Conf. Interval]
-------------+----------------------------------------------------------------
0            |
ipv_version2 |   2.054988   .2680019     5.52   0.000     1.591475      2.6535
ipv_version3 |   .6459268   .0733242    -3.85   0.000     .5170789    .8068817
      female |   1.664704   .1900117     4.47   0.000     1.331005    2.082065
        vid2 |   .8501366   .0082729   -16.68   0.000     .8340757    .8665068
        vid3 |   1.150983   .0059397    27.25   0.000       1.1394    1.162683
        vid4 |   1.032857   .0049758     6.71   0.000      1.02315    1.042655
        vid5 |   .8774016   .0076796   -14.94   0.000     .8624783    .8925831
        vid6 |   1.078792   .0134955     6.06   0.000     1.052663     1.10557
        vid7 |   .8699775   .0097163   -12.47   0.000     .8511408    .8892311
        vid8 |   1.244354   .0101661    26.76   0.000     1.224587     1.26444
       _cons |   1.139277   .0969762     1.53   0.126     .9642159    1.346121
-------------+----------------------------------------------------------------
1            |
ipv_version2 |   2.214611   .3236432     5.44   0.000     1.663041    2.949117
ipv_version3 |    .628279   .0612132    -4.77   0.000     .5190629    .7604752
      female |   1.664704   .1900117     4.47   0.000     1.331005    2.082065
        vid2 |   .8501366   .0082729   -16.68   0.000     .8340757    .8665068
        vid3 |   1.150983   .0059397    27.25   0.000       1.1394    1.162683
        vid4 |   1.032857   .0049758     6.71   0.000      1.02315    1.042655
        vid5 |   .8774016   .0076796   -14.94   0.000     .8624783    .8925831
        vid6 |   1.078792   .0134955     6.06   0.000     1.052663     1.10557
        vid7 |   .8699775   .0097163   -12.47   0.000     .8511408    .8892311
        vid8 |   1.244354   .0101661    26.76   0.000     1.224587     1.26444
       _cons |   .6020262   .0350035    -8.73   0.000     .5371853    .6746937
-------------+----------------------------------------------------------------
2            |
ipv_version2 |   2.846881   .3680378     8.09   0.000     2.209673    3.667841
ipv_version3 |   .5629626   .0725103    -4.46   0.000     .4373646    .7246285
      female |   1.664704   .1900117     4.47   0.000     1.331005    2.082065
        vid2 |   .8501366   .0082729   -16.68   0.000     .8340757    .8665068
        vid3 |   1.150983   .0059397    27.25   0.000       1.1394    1.162683
        vid4 |   1.032857   .0049758     6.71   0.000      1.02315    1.042655
        vid5 |   .8774016   .0076796   -14.94   0.000     .8624783    .8925831
        vid6 |   1.078792   .0134955     6.06   0.000     1.052663     1.10557
        vid7 |   .8699775   .0097163   -12.47   0.000     .8511408    .8892311
        vid8 |   1.244354   .0101661    26.76   0.000     1.224587     1.26444
       _cons |   .3140239   .0191912   -18.95   0.000     .2785752    .3539834
-------------+----------------------------------------------------------------
3            |
ipv_version2 |   3.059949    .321462    10.65   0.000     2.490528    3.759559
ipv_version3 |   .4877618   .0812843    -4.31   0.000     .3518498    .6761739
      female |   1.664704   .1900117     4.47   0.000     1.331005    2.082065
        vid2 |   .8501366   .0082729   -16.68   0.000     .8340757    .8665068
        vid3 |   1.150983   .0059397    27.25   0.000       1.1394    1.162683
        vid4 |   1.032857   .0049758     6.71   0.000      1.02315    1.042655
        vid5 |   .8774016   .0076796   -14.94   0.000     .8624783    .8925831
        vid6 |   1.078792   .0134955     6.06   0.000     1.052663     1.10557
        vid7 |   .8699775   .0097163   -12.47   0.000     .8511408    .8892311
        vid8 |   1.244354   .0101661    26.76   0.000     1.224587     1.26444
       _cons |   .1633242   .0179453   -16.49   0.000     .1316814    .2025706
-------------+----------------------------------------------------------------
4            |
ipv_version2 |   3.513375   .8117898     5.44   0.000     2.233821    5.525869
ipv_version3 |   .5346452   .2257147    -1.48   0.138     .2337266    1.222991
      female |   1.664704   .1900117     4.47   0.000     1.331005    2.082065
        vid2 |   .8501366   .0082729   -16.68   0.000     .8340757    .8665068
        vid3 |   1.150983   .0059397    27.25   0.000       1.1394    1.162683
        vid4 |   1.032857   .0049758     6.71   0.000      1.02315    1.042655
        vid5 |   .8774016   .0076796   -14.94   0.000     .8624783    .8925831
        vid6 |   1.078792   .0134955     6.06   0.000     1.052663     1.10557
        vid7 |   .8699775   .0097163   -12.47   0.000     .8511408    .8892311
        vid8 |   1.244354   .0101661    26.76   0.000     1.224587     1.26444
       _cons |   .0604926   .0128937   -13.16   0.000     .0398358    .0918611
------------------------------------------------------------------------------

. 
. ***** sex stratified estimates: beliefs
. 
. generate femaleXver2=female*ipv_version2

. generate femaleXver3=female*ipv_version3

. 
. * show that interaction terms do not vary across logits
. 
. quietly ologit belief_tot ipv_version2 ipv_version3 female vid2 vid3 vid4 vid5 vid6 vid
> 7 vid8 /*
> */ femaleXver2 femaleXver3, vce(cluster vid)

. brant 

Brant Test of Parallel Regression Assumption

    Variable |      chi2   p>chi2    df
-------------+--------------------------
         All |     78.75    0.003    48
-------------+--------------------------
ipv_version2 |      5.92    0.205     4
ipv_version3 |      6.96    0.138     4
      female |      1.28    0.865     4
        vid2 |      3.31    0.508     4
        vid3 |      3.20    0.525     4
        vid4 |      9.99    0.041     4
        vid5 |      4.45    0.349     4
        vid6 |      2.19    0.701     4
        vid7 |      1.21    0.876     4
        vid8 |      4.78    0.311     4
 femaleXver2 |      3.31    0.507     4
 femaleXver3 |      1.89    0.755     4
----------------------------------------

A significant test statistic provides evidence that the parallel
regression assumption has been violated.

. 
. * beliefs: women
. 
. gologit2 belief_tot ipv_version2 ipv_version3 vid2 vid3 vid4 vid5 vid6 vid7 vid8 if fem
> ale==1, /*
> */ robust cluster(vid) pl(vid2 vid3 vid4 vid5 vid6 vid7 vid8) or

Generalized Ordered Logit Estimates             Number of obs     =        756
                                                Wald chi2(2)      =      17.90
                                                Prob > chi2       =     0.0001
Log pseudolikelihood = -1136.7648               Pseudo R2         =     0.0964

 ( 1)  [0]vid2 - [1]vid2 = 0
 ( 2)  [0]vid3 - [1]vid3 = 0
 ( 3)  [0]vid4 - [1]vid4 = 0
 ( 4)  [0]vid5 - [1]vid5 = 0
 ( 5)  [0]vid6 - [1]vid6 = 0
 ( 6)  [0]vid7 - [1]vid7 = 0
 ( 7)  [0]vid8 - [1]vid8 = 0
 ( 8)  [1]vid2 - [2]vid2 = 0
 ( 9)  [1]vid3 - [2]vid3 = 0
 (10)  [1]vid4 - [2]vid4 = 0
 (11)  [1]vid5 - [2]vid5 = 0
 (12)  [1]vid6 - [2]vid6 = 0
 (13)  [1]vid7 - [2]vid7 = 0
 (14)  [1]vid8 - [2]vid8 = 0
 (15)  [2]vid2 - [3]vid2 = 0
 (16)  [2]vid3 - [3]vid3 = 0
 (17)  [2]vid4 - [3]vid4 = 0
 (18)  [2]vid5 - [3]vid5 = 0
 (19)  [2]vid6 - [3]vid6 = 0
 (20)  [2]vid7 - [3]vid7 = 0
 (21)  [2]vid8 - [3]vid8 = 0
 (22)  [3]vid2 - [4]vid2 = 0
 (23)  [3]vid3 - [4]vid3 = 0
 (24)  [3]vid4 - [4]vid4 = 0
 (25)  [3]vid5 - [4]vid5 = 0
 (26)  [3]vid6 - [4]vid6 = 0
 (27)  [3]vid7 - [4]vid7 = 0
 (28)  [3]vid8 - [4]vid8 = 0
                                    (Std. Err. adjusted for 8 clusters in vid)
------------------------------------------------------------------------------
             |               Robust
  belief_tot | Odds Ratio   Std. Err.      z    P>|z|     [95% Conf. Interval]
-------------+----------------------------------------------------------------
0            |
ipv_version2 |   4.724847    .556128    13.19   0.000     3.751447    5.950817
ipv_version3 |   .5780514   .1367007    -2.32   0.020     .3636384    .9188892
        vid2 |   .8950953   .0082477   -12.03   0.000     .8790752    .9114073
        vid3 |   1.102851   .0103131    10.47   0.000     1.082822    1.123251
        vid4 |    1.03184   .0128217     2.52   0.012     1.007014    1.057279
        vid5 |   .9752608   .0144829    -1.69   0.092     .9472839    1.004064
        vid6 |   .7677749   .0066938   -30.31   0.000     .7547667    .7810073
        vid7 |   .9870154   .0136067    -0.95   0.343     .9607039    1.014048
        vid8 |   1.173754    .019836     9.48   0.000     1.135513    1.213283
       _cons |   1.363797   .1280304     3.31   0.001     1.134595    1.639301
-------------+----------------------------------------------------------------
1            |
ipv_version2 |   5.397789   1.101805     8.26   0.000     3.617989    8.053127
ipv_version3 |    .324847   .0613879    -5.95   0.000     .2242969    .4704729
        vid2 |   .8950953   .0082477   -12.03   0.000     .8790752    .9114073
        vid3 |   1.102851   .0103131    10.47   0.000     1.082822    1.123251
        vid4 |    1.03184   .0128217     2.52   0.012     1.007014    1.057279
        vid5 |   .9752608   .0144829    -1.69   0.092     .9472839    1.004064
        vid6 |   .7677749   .0066938   -30.31   0.000     .7547667    .7810073
        vid7 |   .9870154   .0136067    -0.95   0.343     .9607039    1.014048
        vid8 |   1.173754    .019836     9.48   0.000     1.135513    1.213283
       _cons |   .7599879   .0537174    -3.88   0.000      .661671    .8729136
-------------+----------------------------------------------------------------
2            |
ipv_version2 |   6.315879   1.712282     6.80   0.000     3.712502    10.74486
ipv_version3 |   .3463186   .1026272    -3.58   0.000      .193746    .6190401
        vid2 |   .8950953   .0082477   -12.03   0.000     .8790752    .9114073
        vid3 |   1.102851   .0103131    10.47   0.000     1.082822    1.123251
        vid4 |    1.03184   .0128217     2.52   0.012     1.007014    1.057279
        vid5 |   .9752608   .0144829    -1.69   0.092     .9472839    1.004064
        vid6 |   .7677749   .0066938   -30.31   0.000     .7547667    .7810073
        vid7 |   .9870154   .0136067    -0.95   0.343     .9607039    1.014048
        vid8 |   1.173754    .019836     9.48   0.000     1.135513    1.213283
       _cons |   .3618094   .0549849    -6.69   0.000     .2686093    .4873475
-------------+----------------------------------------------------------------
3            |
ipv_version2 |   4.900495   1.378712     5.65   0.000      2.82333    8.505858
ipv_version3 |   .3394134   .1088169    -3.37   0.001     .1810637    .6362484
        vid2 |   .8950953   .0082477   -12.03   0.000     .8790752    .9114073
        vid3 |   1.102851   .0103131    10.47   0.000     1.082822    1.123251
        vid4 |    1.03184   .0128217     2.52   0.012     1.007014    1.057279
        vid5 |   .9752608   .0144829    -1.69   0.092     .9472839    1.004064
        vid6 |   .7677749   .0066938   -30.31   0.000     .7547667    .7810073
        vid7 |   .9870154   .0136067    -0.95   0.343     .9607039    1.014048
        vid8 |   1.173754    .019836     9.48   0.000     1.135513    1.213283
       _cons |   .1942091   .0357773    -8.90   0.000     .1353508    .2786624
-------------+----------------------------------------------------------------
4            |
ipv_version2 |   3.464202   .7501102     5.74   0.000     2.266162    5.295604
ipv_version3 |   .2202795   .0861941    -3.87   0.000     .1023067    .4742904
        vid2 |   .8950953   .0082477   -12.03   0.000     .8790752    .9114073
        vid3 |   1.102851   .0103131    10.47   0.000     1.082822    1.123251
        vid4 |    1.03184   .0128217     2.52   0.012     1.007014    1.057279
        vid5 |   .9752608   .0144829    -1.69   0.092     .9472839    1.004064
        vid6 |   .7677749   .0066938   -30.31   0.000     .7547667    .7810073
        vid7 |   .9870154   .0136067    -0.95   0.343     .9607039    1.014048
        vid8 |   1.173754    .019836     9.48   0.000     1.135513    1.213283
       _cons |    .095502   .0119265   -18.81   0.000     .0747675    .1219867
------------------------------------------------------------------------------

. 
. * beliefs: men
. 
. gologit2 belief_tot ipv_version2 ipv_version3 vid2 vid3 vid4 vid5 vid6 vid7 vid8 if fem
> ale==0, /*
> */ robust cluster(vid) pl(vid2 vid3 vid4 vid5 vid6 vid7 vid8) or

Generalized Ordered Logit Estimates             Number of obs     =        578
                                                Wald chi2(2)      =      47.05
                                                Prob > chi2       =     0.0000
Log pseudolikelihood = -790.18486               Pseudo R2         =     0.0684

 ( 1)  [0]vid2 - [1]vid2 = 0
 ( 2)  [0]vid3 - [1]vid3 = 0
 ( 3)  [0]vid4 - [1]vid4 = 0
 ( 4)  [0]vid5 - [1]vid5 = 0
 ( 5)  [0]vid6 - [1]vid6 = 0
 ( 6)  [0]vid7 - [1]vid7 = 0
 ( 7)  [0]vid8 - [1]vid8 = 0
 ( 8)  [1]vid2 - [2]vid2 = 0
 ( 9)  [1]vid3 - [2]vid3 = 0
 (10)  [1]vid4 - [2]vid4 = 0
 (11)  [1]vid5 - [2]vid5 = 0
 (12)  [1]vid6 - [2]vid6 = 0
 (13)  [1]vid7 - [2]vid7 = 0
 (14)  [1]vid8 - [2]vid8 = 0
 (15)  [2]vid2 - [3]vid2 = 0
 (16)  [2]vid3 - [3]vid3 = 0
 (17)  [2]vid4 - [3]vid4 = 0
 (18)  [2]vid5 - [3]vid5 = 0
 (19)  [2]vid6 - [3]vid6 = 0
 (20)  [2]vid7 - [3]vid7 = 0
 (21)  [2]vid8 - [3]vid8 = 0
 (22)  [3]vid2 - [4]vid2 = 0
 (23)  [3]vid3 - [4]vid3 = 0
 (24)  [3]vid4 - [4]vid4 = 0
 (25)  [3]vid5 - [4]vid5 = 0
 (26)  [3]vid6 - [4]vid6 = 0
 (27)  [3]vid7 - [4]vid7 = 0
 (28)  [3]vid8 - [4]vid8 = 0
                                    (Std. Err. adjusted for 8 clusters in vid)
------------------------------------------------------------------------------
             |               Robust
  belief_tot | Odds Ratio   Std. Err.      z    P>|z|     [95% Conf. Interval]
-------------+----------------------------------------------------------------
0            |
ipv_version2 |   3.387579   .8947665     4.62   0.000     2.018646    5.684844
ipv_version3 |   .9136059   .1645444    -0.50   0.616     .6418799    1.300361
        vid2 |   .8702666   .0149173    -8.11   0.000     .8415148    .9000007
        vid3 |   1.400896   .0281862    16.75   0.000     1.346727    1.457243
        vid4 |   1.357123   .0304949    13.59   0.000     1.298652    1.418228
        vid5 |    .762462   .0158521   -13.04   0.000      .732017    .7941733
        vid6 |   .8012606   .0140448   -12.64   0.000     .7742007    .8292663
        vid7 |   1.204706   .0176009    12.75   0.000     1.170698    1.239702
        vid8 |   1.071504    .008503     8.70   0.000     1.054967      1.0883
       _cons |   .6529917   .0845377    -3.29   0.001     .5066515    .8416005
-------------+----------------------------------------------------------------
1            |
ipv_version2 |   4.100018   .9491599     6.09   0.000     2.604543    6.454166
ipv_version3 |   .5983597   .1537133    -2.00   0.046     .3616565    .9899845
        vid2 |   .8702666   .0149173    -8.11   0.000     .8415148    .9000007
        vid3 |   1.400896   .0281862    16.75   0.000     1.346727    1.457243
        vid4 |   1.357123   .0304949    13.59   0.000     1.298652    1.418228
        vid5 |    .762462   .0158521   -13.04   0.000      .732017    .7941733
        vid6 |   .8012606   .0140448   -12.64   0.000     .7742007    .8292663
        vid7 |   1.204706   .0176009    12.75   0.000     1.170698    1.239702
        vid8 |   1.071504    .008503     8.70   0.000     1.054967      1.0883
       _cons |   .3516766   .0453622    -8.10   0.000     .2731168    .4528335
-------------+----------------------------------------------------------------
2            |
ipv_version2 |   5.578659   1.312906     7.30   0.000      3.51726    8.848207
ipv_version3 |   .6859868   .2144943    -1.21   0.228      .371674    1.266104
        vid2 |   .8702666   .0149173    -8.11   0.000     .8415148    .9000007
        vid3 |   1.400896   .0281862    16.75   0.000     1.346727    1.457243
        vid4 |   1.357123   .0304949    13.59   0.000     1.298652    1.418228
        vid5 |    .762462   .0158521   -13.04   0.000      .732017    .7941733
        vid6 |   .8012606   .0140448   -12.64   0.000     .7742007    .8292663
        vid7 |   1.204706   .0176009    12.75   0.000     1.170698    1.239702
        vid8 |   1.071504    .008503     8.70   0.000     1.054967      1.0883
       _cons |   .1590128   .0351834    -8.31   0.000      .103061    .2453409
-------------+----------------------------------------------------------------
3            |
ipv_version2 |   6.119982   1.273585     8.71   0.000     4.070173    9.202112
ipv_version3 |   .4083937   .1753862    -2.09   0.037     .1760059    .9476125
        vid2 |   .8702666   .0149173    -8.11   0.000     .8415148    .9000007
        vid3 |   1.400896   .0281862    16.75   0.000     1.346727    1.457243
        vid4 |   1.357123   .0304949    13.59   0.000     1.298652    1.418228
        vid5 |    .762462   .0158521   -13.04   0.000      .732017    .7941733
        vid6 |   .8012606   .0140448   -12.64   0.000     .7742007    .8292663
        vid7 |   1.204706   .0176009    12.75   0.000     1.170698    1.239702
        vid8 |   1.071504    .008503     8.70   0.000     1.054967      1.0883
       _cons |   .0864419   .0205352   -10.31   0.000      .054264    .1377008
-------------+----------------------------------------------------------------
4            |
ipv_version2 |   5.751735   3.475196     2.90   0.004     1.759962    18.79726
ipv_version3 |   .3412529   .4271674    -0.86   0.390     .0293477    3.968059
        vid2 |   .8702666   .0149173    -8.11   0.000     .8415148    .9000007
        vid3 |   1.400896   .0281862    16.75   0.000     1.346727    1.457243
        vid4 |   1.357123   .0304949    13.59   0.000     1.298652    1.418228
        vid5 |    .762462   .0158521   -13.04   0.000      .732017    .7941733
        vid6 |   .8012606   .0140448   -12.64   0.000     .7742007    .8292663
        vid7 |   1.204706   .0176009    12.75   0.000     1.170698    1.239702
        vid8 |   1.071504    .008503     8.70   0.000     1.054967      1.0883
       _cons |   .0287219   .0158992    -6.41   0.000     .0097056    .0849967
------------------------------------------------------------------------------

. 
. * test interaction term
. 
. xi: gologit2 belief_tot ipv_version2 ipv_version3 female vid2 vid3 vid4 vid5 vid6 vid7 
> vid8 /*
> */ femaleXver2 femaleXver3, pl(female vid2 vid3 vid4 vid5 vid6 vid7 vid8 femaleXver2 fe
> maleXver3) /*
> */ robust cluster(vid) or 

Generalized Ordered Logit Estimates             Number of obs     =      1,334
                                                Wald chi2(2)      =       5.55
                                                Prob > chi2       =     0.0625
Log pseudolikelihood = -1933.2965               Pseudo R2         =     0.0888

 ( 1)  [0]female - [1]female = 0
 ( 2)  [0]vid2 - [1]vid2 = 0
 ( 3)  [0]vid3 - [1]vid3 = 0
 ( 4)  [0]vid4 - [1]vid4 = 0
 ( 5)  [0]vid5 - [1]vid5 = 0
 ( 6)  [0]vid6 - [1]vid6 = 0
 ( 7)  [0]vid7 - [1]vid7 = 0
 ( 8)  [0]vid8 - [1]vid8 = 0
 ( 9)  [0]femaleXver2 - [1]femaleXver2 = 0
 (10)  [0]femaleXver3 - [1]femaleXver3 = 0
 (11)  [1]female - [2]female = 0
 (12)  [1]vid2 - [2]vid2 = 0
 (13)  [1]vid3 - [2]vid3 = 0
 (14)  [1]vid4 - [2]vid4 = 0
 (15)  [1]vid5 - [2]vid5 = 0
 (16)  [1]vid6 - [2]vid6 = 0
 (17)  [1]vid7 - [2]vid7 = 0
 (18)  [1]vid8 - [2]vid8 = 0
 (19)  [1]femaleXver2 - [2]femaleXver2 = 0
 (20)  [1]femaleXver3 - [2]femaleXver3 = 0
 (21)  [2]female - [3]female = 0
 (22)  [2]vid2 - [3]vid2 = 0
 (23)  [2]vid3 - [3]vid3 = 0
 (24)  [2]vid4 - [3]vid4 = 0
 (25)  [2]vid5 - [3]vid5 = 0
 (26)  [2]vid6 - [3]vid6 = 0
 (27)  [2]vid7 - [3]vid7 = 0
 (28)  [2]vid8 - [3]vid8 = 0
 (29)  [2]femaleXver2 - [3]femaleXver2 = 0
 (30)  [2]femaleXver3 - [3]femaleXver3 = 0
 (31)  [3]female - [4]female = 0
 (32)  [3]vid2 - [4]vid2 = 0
 (33)  [3]vid3 - [4]vid3 = 0
 (34)  [3]vid4 - [4]vid4 = 0
 (35)  [3]vid5 - [4]vid5 = 0
 (36)  [3]vid6 - [4]vid6 = 0
 (37)  [3]vid7 - [4]vid7 = 0
 (38)  [3]vid8 - [4]vid8 = 0
 (39)  [3]femaleXver2 - [4]femaleXver2 = 0
 (40)  [3]femaleXver3 - [4]femaleXver3 = 0
                                    (Std. Err. adjusted for 8 clusters in vid)
------------------------------------------------------------------------------
             |               Robust
  belief_tot | Odds Ratio   Std. Err.      z    P>|z|     [95% Conf. Interval]
-------------+----------------------------------------------------------------
0            |
ipv_version2 |   3.754475    1.04393     4.76   0.000     2.177073    6.474787
ipv_version3 |   .9314182   .1736636    -0.38   0.703      .646305    1.342307
      female |   2.022897   .2997926     4.75   0.000     1.512955    2.704715
        vid2 |   .9001344   .0094121   -10.06   0.000     .8818747    .9187722
        vid3 |   1.233414   .0111655    23.17   0.000     1.211723    1.255493
        vid4 |    1.17349   .0078797    23.83   0.000     1.158147    1.189036
        vid5 |   .8940951   .0091197   -10.97   0.000     .8763983    .9121493
        vid6 |   .7847265   .0108791   -17.49   0.000      .763691    .8063414
        vid7 |   1.091396   .0060637    15.74   0.000     1.079575    1.103345
        vid8 |   1.141737   .0091649    16.51   0.000     1.123914    1.159842
 femaleXver2 |   1.074208   .3619602     0.21   0.832     .5549726    2.079244
 femaleXver3 |   .6063882   .0785425    -3.86   0.000     .4704341    .7816326
       _cons |   .6625135   .0849119    -3.21   0.001     .5153465    .8517069
-------------+----------------------------------------------------------------
1            |
ipv_version2 |   4.417715   .9930541     6.61   0.000     2.843524    6.863388
ipv_version3 |   .5580145   .1014712    -3.21   0.001     .3907151    .7969495
      female |   2.022897   .2997926     4.75   0.000     1.512955    2.704715
        vid2 |   .9001344   .0094121   -10.06   0.000     .8818747    .9187722
        vid3 |   1.233414   .0111655    23.17   0.000     1.211723    1.255493
        vid4 |    1.17349   .0078797    23.83   0.000     1.158147    1.189036
        vid5 |   .8940951   .0091197   -10.97   0.000     .8763983    .9121493
        vid6 |   .7847265   .0108791   -17.49   0.000      .763691    .8063414
        vid7 |   1.091396   .0060637    15.74   0.000     1.079575    1.103345
        vid8 |   1.141737   .0091649    16.51   0.000     1.123914    1.159842
 femaleXver2 |   1.074208   .3619602     0.21   0.832     .5549726    2.079244
 femaleXver3 |   .6063882   .0785425    -3.86   0.000     .4704341    .7816326
       _cons |    .364677   .0393209    -9.36   0.000     .2952083    .4504932
-------------+----------------------------------------------------------------
2            |
ipv_version2 |   5.539654   .7410952    12.80   0.000     4.261955    7.200397
ipv_version3 |   .6102884   .1550424    -1.94   0.052     .3709277    1.004109
      female |   2.022897   .2997926     4.75   0.000     1.512955    2.704715
        vid2 |   .9001344   .0094121   -10.06   0.000     .8818747    .9187722
        vid3 |   1.233414   .0111655    23.17   0.000     1.211723    1.255493
        vid4 |    1.17349   .0078797    23.83   0.000     1.158147    1.189036
        vid5 |   .8940951   .0091197   -10.97   0.000     .8763983    .9121493
        vid6 |   .7847265   .0108791   -17.49   0.000      .763691    .8063414
        vid7 |   1.091396   .0060637    15.74   0.000     1.079575    1.103345
        vid8 |   1.141737   .0091649    16.51   0.000     1.123914    1.159842
 femaleXver2 |   1.074208   .3619602     0.21   0.832     .5549726    2.079244
 femaleXver3 |   .6063882   .0785425    -3.86   0.000     .4704341    .7816326
       _cons |   .1709664   .0183944   -16.42   0.000     .1384616     .211102
-------------+----------------------------------------------------------------
3            |
ipv_version2 |   4.980343   .8081914     9.89   0.000      3.62351    6.845246
ipv_version3 |   .4957777   .1566815    -2.22   0.026      .266858    .9210723
      female |   2.022897   .2997926     4.75   0.000     1.512955    2.704715
        vid2 |   .9001344   .0094121   -10.06   0.000     .8818747    .9187722
        vid3 |   1.233414   .0111655    23.17   0.000     1.211723    1.255493
        vid4 |    1.17349   .0078797    23.83   0.000     1.158147    1.189036
        vid5 |   .8940951   .0091197   -10.97   0.000     .8763983    .9121493
        vid6 |   .7847265   .0108791   -17.49   0.000      .763691    .8063414
        vid7 |   1.091396   .0060637    15.74   0.000     1.079575    1.103345
        vid8 |   1.141737   .0091649    16.51   0.000     1.123914    1.159842
 femaleXver2 |   1.074208   .3619602     0.21   0.832     .5549726    2.079244
 femaleXver3 |   .6063882   .0785425    -3.86   0.000     .4704341    .7816326
       _cons |   .0923785   .0108279   -20.32   0.000     .0734175    .1162363
-------------+----------------------------------------------------------------
4            |
ipv_version2 |   3.734817   1.209479     4.07   0.000     1.979781    7.045656
ipv_version3 |    .341546   .1455987    -2.52   0.012     .1481112    .7876089
      female |   2.022897   .2997926     4.75   0.000     1.512955    2.704715
        vid2 |   .9001344   .0094121   -10.06   0.000     .8818747    .9187722
        vid3 |   1.233414   .0111655    23.17   0.000     1.211723    1.255493
        vid4 |    1.17349   .0078797    23.83   0.000     1.158147    1.189036
        vid5 |   .8940951   .0091197   -10.97   0.000     .8763983    .9121493
        vid6 |   .7847265   .0108791   -17.49   0.000      .763691    .8063414
        vid7 |   1.091396   .0060637    15.74   0.000     1.079575    1.103345
        vid8 |   1.141737   .0091649    16.51   0.000     1.123914    1.159842
 femaleXver2 |   1.074208   .3619602     0.21   0.832     .5549726    2.079244
 femaleXver3 |   .6063882   .0785425    -3.86   0.000     .4704341    .7816326
       _cons |   .0410793   .0078645   -16.67   0.000     .0282269    .0597836
------------------------------------------------------------------------------

. test femaleXver2 femaleXver3

 ( 1)  [0]femaleXver2 = 0
 ( 2)  [1]femaleXver2 = 0
 ( 3)  [2]femaleXver2 = 0
 ( 4)  [3]femaleXver2 = 0
 ( 5)  [4]femaleXver2 = 0
 ( 6)  [0]femaleXver3 = 0
 ( 7)  [1]femaleXver3 = 0
 ( 8)  [2]femaleXver3 = 0
 ( 9)  [3]femaleXver3 = 0
 (10)  [4]femaleXver3 = 0
       Constraint 1 dropped
       Constraint 2 dropped
       Constraint 3 dropped
       Constraint 5 dropped
       Constraint 7 dropped
       Constraint 8 dropped
       Constraint 9 dropped
       Constraint 10 dropped

           chi2(  2) =   14.93
         Prob > chi2 =    0.0006

. 
. ***** sex stratified estimates: norms
. 
. * show that interaction terms do not vary across logits
. 
. quietly ologit norm_tot ipv_version2 ipv_version3 female vid2 vid3 vid4 vid5 vid6 vid7 
> vid8 /*
> */ femaleXver2 femaleXver3, vce(cluster vid)

. brant 

Brant Test of Parallel Regression Assumption

    Variable |      chi2   p>chi2    df
-------------+--------------------------
         All |     44.17    0.630    48
-------------+--------------------------
ipv_version2 |      7.61    0.107     4
ipv_version3 |      0.72    0.949     4
      female |      1.11    0.893     4
        vid2 |      2.26    0.689     4
        vid3 |      2.26    0.688     4
        vid4 |      0.20    0.995     4
        vid5 |      1.76    0.780     4
        vid6 |      1.97    0.740     4
        vid7 |      0.51    0.972     4
        vid8 |      8.61    0.072     4
 femaleXver2 |      1.92    0.751     4
 femaleXver3 |      3.46    0.484     4
----------------------------------------

A significant test statistic provides evidence that the parallel
regression assumption has been violated.

. 
. * norms: women
. 
. gologit2 norm_tot ipv_version2 ipv_version3 vid2 vid3 vid4 vid5 vid6 vid7 vid8 if femal
> e==1, /*
> */ robust cluster(vid) pl(vid2 vid3 vid4 vid5 vid6 vid7 vid8) or

Generalized Ordered Logit Estimates             Number of obs     =        756
                                                Wald chi2(2)      =      11.42
                                                Prob > chi2       =     0.0033
Log pseudolikelihood = -1229.3564               Pseudo R2         =     0.0509

 ( 1)  [0]vid2 - [1]vid2 = 0
 ( 2)  [0]vid3 - [1]vid3 = 0
 ( 3)  [0]vid4 - [1]vid4 = 0
 ( 4)  [0]vid5 - [1]vid5 = 0
 ( 5)  [0]vid6 - [1]vid6 = 0
 ( 6)  [0]vid7 - [1]vid7 = 0
 ( 7)  [0]vid8 - [1]vid8 = 0
 ( 8)  [1]vid2 - [2]vid2 = 0
 ( 9)  [1]vid3 - [2]vid3 = 0
 (10)  [1]vid4 - [2]vid4 = 0
 (11)  [1]vid5 - [2]vid5 = 0
 (12)  [1]vid6 - [2]vid6 = 0
 (13)  [1]vid7 - [2]vid7 = 0
 (14)  [1]vid8 - [2]vid8 = 0
 (15)  [2]vid2 - [3]vid2 = 0
 (16)  [2]vid3 - [3]vid3 = 0
 (17)  [2]vid4 - [3]vid4 = 0
 (18)  [2]vid5 - [3]vid5 = 0
 (19)  [2]vid6 - [3]vid6 = 0
 (20)  [2]vid7 - [3]vid7 = 0
 (21)  [2]vid8 - [3]vid8 = 0
 (22)  [3]vid2 - [4]vid2 = 0
 (23)  [3]vid3 - [4]vid3 = 0
 (24)  [3]vid4 - [4]vid4 = 0
 (25)  [3]vid5 - [4]vid5 = 0
 (26)  [3]vid6 - [4]vid6 = 0
 (27)  [3]vid7 - [4]vid7 = 0
 (28)  [3]vid8 - [4]vid8 = 0
                                    (Std. Err. adjusted for 8 clusters in vid)
------------------------------------------------------------------------------
             |               Robust
    norm_tot | Odds Ratio   Std. Err.      z    P>|z|     [95% Conf. Interval]
-------------+----------------------------------------------------------------
0            |
ipv_version2 |   3.162475   .3888262     9.36   0.000      2.48526    4.024225
ipv_version3 |    .750477   .1130768    -1.91   0.057     .5585784    1.008302
        vid2 |   .8955497   .0093613   -10.55   0.000     .8773886    .9140868
        vid3 |   .8576549   .0136776    -9.63   0.000     .8312619     .884886
        vid4 |   .9193266   .0124133    -6.23   0.000     .8953162     .943981
        vid5 |   .9294111   .0187233    -3.63   0.000     .8934291    .9668422
        vid6 |   .9137324   .0075719   -10.89   0.000     .8990116    .9286943
        vid7 |   .6809129   .0139434   -18.77   0.000     .6541255    .7087973
        vid8 |   1.032857   .0149946     2.23   0.026     1.003883    1.062668
       _cons |   1.806142   .1231226     8.67   0.000     1.580252    2.064321
-------------+----------------------------------------------------------------
1            |
ipv_version2 |   3.056565   .7750273     4.41   0.000     1.859523    5.024185
ipv_version3 |   .7104015   .1138632    -2.13   0.033     .5188875    .9726006
        vid2 |   .8955497   .0093613   -10.55   0.000     .8773886    .9140868
        vid3 |   .8576549   .0136776    -9.63   0.000     .8312619     .884886
        vid4 |   .9193266   .0124133    -6.23   0.000     .8953162     .943981
        vid5 |   .9294111   .0187233    -3.63   0.000     .8934291    .9668422
        vid6 |   .9137324   .0075719   -10.89   0.000     .8990116    .9286943
        vid7 |   .6809129   .0139434   -18.77   0.000     .6541255    .7087973
        vid8 |   1.032857   .0149946     2.23   0.026     1.003883    1.062668
       _cons |   1.001997   .1410875     0.01   0.989     .7603472    1.320447
-------------+----------------------------------------------------------------
2            |
ipv_version2 |   3.721564   .7009388     6.98   0.000     2.572798    5.383259
ipv_version3 |   .5418863   .1089025    -3.05   0.002     .3654624    .8034772
        vid2 |   .8955497   .0093613   -10.55   0.000     .8773886    .9140868
        vid3 |   .8576549   .0136776    -9.63   0.000     .8312619     .884886
        vid4 |   .9193266   .0124133    -6.23   0.000     .8953162     .943981
        vid5 |   .9294111   .0187233    -3.63   0.000     .8934291    .9668422
        vid6 |   .9137324   .0075719   -10.89   0.000     .8990116    .9286943
        vid7 |   .6809129   .0139434   -18.77   0.000     .6541255    .7087973
        vid8 |   1.032857   .0149946     2.23   0.026     1.003883    1.062668
       _cons |   .5273799   .0566891    -5.95   0.000     .4271951    .6510598
-------------+----------------------------------------------------------------
3            |
ipv_version2 |    3.68926   .6114989     7.88   0.000     2.665948    5.105366
ipv_version3 |   .4501888   .0956371    -3.76   0.000     .2968707    .6826877
        vid2 |   .8955497   .0093613   -10.55   0.000     .8773886    .9140868
        vid3 |   .8576549   .0136776    -9.63   0.000     .8312619     .884886
        vid4 |   .9193266   .0124133    -6.23   0.000     .8953162     .943981
        vid5 |   .9294111   .0187233    -3.63   0.000     .8934291    .9668422
        vid6 |   .9137324   .0075719   -10.89   0.000     .8990116    .9286943
        vid7 |   .6809129   .0139434   -18.77   0.000     .6541255    .7087973
        vid8 |   1.032857   .0149946     2.23   0.026     1.003883    1.062668
       _cons |   .2869266    .027312   -13.12   0.000     .2380929    .3457762
-------------+----------------------------------------------------------------
4            |
ipv_version2 |   4.097046    1.00967     5.72   0.000     2.527564    6.641094
ipv_version3 |   .4954439   .2144688    -1.62   0.105     .2120921    1.157349
        vid2 |   .8955497   .0093613   -10.55   0.000     .8773886    .9140868
        vid3 |   .8576549   .0136776    -9.63   0.000     .8312619     .884886
        vid4 |   .9193266   .0124133    -6.23   0.000     .8953162     .943981
        vid5 |   .9294111   .0187233    -3.63   0.000     .8934291    .9668422
        vid6 |   .9137324   .0075719   -10.89   0.000     .8990116    .9286943
        vid7 |   .6809129   .0139434   -18.77   0.000     .6541255    .7087973
        vid8 |   1.032857   .0149946     2.23   0.026     1.003883    1.062668
       _cons |   .1032702   .0177753   -13.19   0.000     .0736992    .1447063
------------------------------------------------------------------------------

. 
. * norms: men
. 
. gologit2 norm_tot ipv_version2 ipv_version3 vid2 vid3 vid4 vid5 vid6 vid7 vid8 if femal
> e==0, /*
> */ robust cluster(vid) pl(vid2 vid3 vid4 vid5 vid6 vid7 vid8) or /* men */

Generalized Ordered Logit Estimates             Number of obs     =        578
                                                Wald chi2(2)      =       1.49
                                                Prob > chi2       =     0.4752
Log pseudolikelihood = -880.70124               Pseudo R2         =     0.0293

 ( 1)  [0]vid2 - [1]vid2 = 0
 ( 2)  [0]vid3 - [1]vid3 = 0
 ( 3)  [0]vid4 - [1]vid4 = 0
 ( 4)  [0]vid5 - [1]vid5 = 0
 ( 5)  [0]vid6 - [1]vid6 = 0
 ( 6)  [0]vid7 - [1]vid7 = 0
 ( 7)  [0]vid8 - [1]vid8 = 0
 ( 8)  [1]vid2 - [2]vid2 = 0
 ( 9)  [1]vid3 - [2]vid3 = 0
 (10)  [1]vid4 - [2]vid4 = 0
 (11)  [1]vid5 - [2]vid5 = 0
 (12)  [1]vid6 - [2]vid6 = 0
 (13)  [1]vid7 - [2]vid7 = 0
 (14)  [1]vid8 - [2]vid8 = 0
 (15)  [2]vid2 - [3]vid2 = 0
 (16)  [2]vid3 - [3]vid3 = 0
 (17)  [2]vid4 - [3]vid4 = 0
 (18)  [2]vid5 - [3]vid5 = 0
 (19)  [2]vid6 - [3]vid6 = 0
 (20)  [2]vid7 - [3]vid7 = 0
 (21)  [2]vid8 - [3]vid8 = 0
 (22)  [3]vid2 - [4]vid2 = 0
 (23)  [3]vid3 - [4]vid3 = 0
 (24)  [3]vid4 - [4]vid4 = 0
 (25)  [3]vid5 - [4]vid5 = 0
 (26)  [3]vid6 - [4]vid6 = 0
 (27)  [3]vid7 - [4]vid7 = 0
 (28)  [3]vid8 - [4]vid8 = 0
                                    (Std. Err. adjusted for 8 clusters in vid)
------------------------------------------------------------------------------
             |               Robust
    norm_tot | Odds Ratio   Std. Err.      z    P>|z|     [95% Conf. Interval]
-------------+----------------------------------------------------------------
0            |
ipv_version2 |   1.307392   .4092899     0.86   0.392     .7078321    2.414802
ipv_version3 |   .5387391   .0874513    -3.81   0.000     .3919282    .7405434
        vid2 |   .7346009   .0142558   -15.89   0.000     .7071847      .76308
        vid3 |   1.661149   .0510676    16.51   0.000     1.564014    1.764317
        vid4 |    1.18371   .0205318     9.72   0.000     1.144145    1.224644
        vid5 |   .8085436   .0142411   -12.07   0.000     .7811078     .836943
        vid6 |   1.346244   .0173882    23.02   0.000     1.312592     1.38076
        vid7 |   1.192581   .0275769     7.62   0.000     1.139738    1.247874
        vid8 |   1.587571   .0470416    15.60   0.000     1.497998    1.682501
       _cons |   1.215473   .2043913     1.16   0.246     .8741962    1.689979
-------------+----------------------------------------------------------------
1            |
ipv_version2 |   1.526924   .3498151     1.85   0.065      .974561    2.392357
ipv_version3 |    .543815    .122532    -2.70   0.007     .3496705    .8457526
        vid2 |   .7346009   .0142558   -15.89   0.000     .7071847      .76308
        vid3 |   1.661149   .0510676    16.51   0.000     1.564014    1.764317
        vid4 |    1.18371   .0205318     9.72   0.000     1.144145    1.224644
        vid5 |   .8085436   .0142411   -12.07   0.000     .7811078     .836943
        vid6 |   1.346244   .0173882    23.02   0.000     1.312592     1.38076
        vid7 |   1.192581   .0275769     7.62   0.000     1.139738    1.247874
        vid8 |   1.587571   .0470416    15.60   0.000     1.497998    1.682501
       _cons |   .6031025    .069069    -4.42   0.000     .4818469    .7548718
-------------+----------------------------------------------------------------
2            |
ipv_version2 |   2.045462   .3272339     4.47   0.000     1.494913    2.798769
ipv_version3 |   .6237061   .1617566    -1.82   0.069     .3751657      1.0369
        vid2 |   .7346009   .0142558   -15.89   0.000     .7071847      .76308
        vid3 |   1.661149   .0510676    16.51   0.000     1.564014    1.764317
        vid4 |    1.18371   .0205318     9.72   0.000     1.144145    1.224644
        vid5 |   .8085436   .0142411   -12.07   0.000     .7811078     .836943
        vid6 |   1.346244   .0173882    23.02   0.000     1.312592     1.38076
        vid7 |   1.192581   .0275769     7.62   0.000     1.139738    1.247874
        vid8 |   1.587571   .0470416    15.60   0.000     1.497998    1.682501
       _cons |   .3099187   .0432993    -8.38   0.000     .2356812    .4075403
-------------+----------------------------------------------------------------
3            |
ipv_version2 |   2.399969   .6097059     3.45   0.001      1.45868    3.948673
ipv_version3 |   .5773202   .2725396    -1.16   0.245     .2288675    1.456295
        vid2 |   .7346009   .0142558   -15.89   0.000     .7071847      .76308
        vid3 |   1.661149   .0510676    16.51   0.000     1.564014    1.764317
        vid4 |    1.18371   .0205318     9.72   0.000     1.144145    1.224644
        vid5 |   .8085436   .0142411   -12.07   0.000     .7811078     .836943
        vid6 |   1.346244   .0173882    23.02   0.000     1.312592     1.38076
        vid7 |   1.192581   .0275769     7.62   0.000     1.139738    1.247874
        vid8 |   1.587571   .0470416    15.60   0.000     1.497998    1.682501
       _cons |   .1496296   .0347146    -8.19   0.000     .0949589    .2357757
-------------+----------------------------------------------------------------
4            |
ipv_version2 |   2.849246   1.059113     2.82   0.005     1.375067     5.90386
ipv_version3 |    .624183   .4570893    -0.64   0.520     .1485857    2.622086
        vid2 |   .7346009   .0142558   -15.89   0.000     .7071847      .76308
        vid3 |   1.661149   .0510676    16.51   0.000     1.564014    1.764317
        vid4 |    1.18371   .0205318     9.72   0.000     1.144145    1.224644
        vid5 |   .8085436   .0142411   -12.07   0.000     .7811078     .836943
        vid6 |   1.346244   .0173882    23.02   0.000     1.312592     1.38076
        vid7 |   1.192581   .0275769     7.62   0.000     1.139738    1.247874
        vid8 |   1.587571   .0470416    15.60   0.000     1.497998    1.682501
       _cons |   .0580109   .0276696    -5.97   0.000     .0227777    .1477439
------------------------------------------------------------------------------

. 
. * test interaction term
. 
. xi: gologit2 norm_tot ipv_version2 ipv_version3 female vid2 vid3 vid4 vid5 vid6 vid7 vi
> d8 /*
> */ femaleXver2 femaleXver3, pl(female vid2 vid3 vid4 vid5 vid6 vid7 vid8 femaleXver2 fe
> maleXver3) /*
> */ robust cluster(vid) or 

Generalized Ordered Logit Estimates             Number of obs     =      1,334
                                                Wald chi2(2)      =       5.17
                                                Prob > chi2       =     0.0752
Log pseudolikelihood =  -2118.116               Pseudo R2         =     0.0443

 ( 1)  [0]female - [1]female = 0
 ( 2)  [0]vid2 - [1]vid2 = 0
 ( 3)  [0]vid3 - [1]vid3 = 0
 ( 4)  [0]vid4 - [1]vid4 = 0
 ( 5)  [0]vid5 - [1]vid5 = 0
 ( 6)  [0]vid6 - [1]vid6 = 0
 ( 7)  [0]vid7 - [1]vid7 = 0
 ( 8)  [0]vid8 - [1]vid8 = 0
 ( 9)  [0]femaleXver2 - [1]femaleXver2 = 0
 (10)  [0]femaleXver3 - [1]femaleXver3 = 0
 (11)  [1]female - [2]female = 0
 (12)  [1]vid2 - [2]vid2 = 0
 (13)  [1]vid3 - [2]vid3 = 0
 (14)  [1]vid4 - [2]vid4 = 0
 (15)  [1]vid5 - [2]vid5 = 0
 (16)  [1]vid6 - [2]vid6 = 0
 (17)  [1]vid7 - [2]vid7 = 0
 (18)  [1]vid8 - [2]vid8 = 0
 (19)  [1]femaleXver2 - [2]femaleXver2 = 0
 (20)  [1]femaleXver3 - [2]femaleXver3 = 0
 (21)  [2]female - [3]female = 0
 (22)  [2]vid2 - [3]vid2 = 0
 (23)  [2]vid3 - [3]vid3 = 0
 (24)  [2]vid4 - [3]vid4 = 0
 (25)  [2]vid5 - [3]vid5 = 0
 (26)  [2]vid6 - [3]vid6 = 0
 (27)  [2]vid7 - [3]vid7 = 0
 (28)  [2]vid8 - [3]vid8 = 0
 (29)  [2]femaleXver2 - [3]femaleXver2 = 0
 (30)  [2]femaleXver3 - [3]femaleXver3 = 0
 (31)  [3]female - [4]female = 0
 (32)  [3]vid2 - [4]vid2 = 0
 (33)  [3]vid3 - [4]vid3 = 0
 (34)  [3]vid4 - [4]vid4 = 0
 (35)  [3]vid5 - [4]vid5 = 0
 (36)  [3]vid6 - [4]vid6 = 0
 (37)  [3]vid7 - [4]vid7 = 0
 (38)  [3]vid8 - [4]vid8 = 0
 (39)  [3]femaleXver2 - [4]femaleXver2 = 0
 (40)  [3]femaleXver3 - [4]femaleXver3 = 0
                                    (Std. Err. adjusted for 8 clusters in vid)
------------------------------------------------------------------------------
             |               Robust
    norm_tot | Odds Ratio   Std. Err.      z    P>|z|     [95% Conf. Interval]
-------------+----------------------------------------------------------------
0            |
ipv_version2 |   1.451475   .3909311     1.38   0.167     .8561548    2.460746
ipv_version3 |   .5807836    .084213    -3.75   0.000     .4371103    .7716807
      female |   1.254867   .2608198     1.09   0.275     .8349842    1.885893
        vid2 |   .8486625   .0086541   -16.09   0.000     .8318693    .8657948
        vid3 |   1.160412   .0090696    19.04   0.000     1.142771    1.178325
        vid4 |   1.038661   .0065391     6.03   0.000     1.025924    1.051557
        vid5 |   .8948181   .0124975    -7.96   0.000     .8706556    .9196512
        vid6 |   1.080876   .0134682     6.24   0.000     1.054798    1.107598
        vid7 |   .8808432   .0120414    -9.28   0.000     .8575557     .904763
        vid8 |   1.257772   .0122384    23.57   0.000     1.234012    1.281989
 femaleXver2 |   1.919249   .6177648     2.03   0.043     1.021288    3.606739
 femaleXver3 |   1.217535   .3518236     0.68   0.496     .6910587    2.145102
       _cons |    1.32325   .1823934     2.03   0.042     1.009984    1.733682
-------------+----------------------------------------------------------------
1            |
ipv_version2 |   1.539173   .3145999     2.11   0.035     1.031112    2.297572
ipv_version3 |   .5643722   .1274179    -2.53   0.011     .3625688    .8784978
      female |   1.254867   .2608198     1.09   0.275     .8349842    1.885893
        vid2 |   .8486625   .0086541   -16.09   0.000     .8318693    .8657948
        vid3 |   1.160412   .0090696    19.04   0.000     1.142771    1.178325
        vid4 |   1.038661   .0065391     6.03   0.000     1.025924    1.051557
        vid5 |   .8948181   .0124975    -7.96   0.000     .8706556    .9196512
        vid6 |   1.080876   .0134682     6.24   0.000     1.054798    1.107598
        vid7 |   .8808432   .0120414    -9.28   0.000     .8575557     .904763
        vid8 |   1.257772   .0122384    23.57   0.000     1.234012    1.281989
 femaleXver2 |   1.919249   .6177648     2.03   0.043     1.021288    3.606739
 femaleXver3 |   1.217535   .3518236     0.68   0.496     .6910587    2.145102
       _cons |   .7015422   .0736368    -3.38   0.001     .5710946    .8617861
-------------+----------------------------------------------------------------
2            |
ipv_version2 |   1.952594   .3586469     3.64   0.000     1.362278    2.798712
ipv_version3 |   .5045795   .1188119    -2.90   0.004      .318053    .8004969
      female |   1.254867   .2608198     1.09   0.275     .8349842    1.885893
        vid2 |   .8486625   .0086541   -16.09   0.000     .8318693    .8657948
        vid3 |   1.160412   .0090696    19.04   0.000     1.142771    1.178325
        vid4 |   1.038661   .0065391     6.03   0.000     1.025924    1.051557
        vid5 |   .8948181   .0124975    -7.96   0.000     .8706556    .9196512
        vid6 |   1.080876   .0134682     6.24   0.000     1.054798    1.107598
        vid7 |   .8808432   .0120414    -9.28   0.000     .8575557     .904763
        vid8 |   1.257772   .0122384    23.57   0.000     1.234012    1.281989
 femaleXver2 |   1.919249   .6177648     2.03   0.043     1.021288    3.606739
 femaleXver3 |   1.217535   .3518236     0.68   0.496     .6910587    2.145102
       _cons |   .3673823   .0394686    -9.32   0.000     .2976268    .4534866
-------------+----------------------------------------------------------------
3            |
ipv_version2 |   2.065068   .4124043     3.63   0.000     1.396192    3.054383
ipv_version3 |   .4358276   .1290015    -2.81   0.005     .2439866     .778509
      female |   1.254867   .2608198     1.09   0.275     .8349842    1.885893
        vid2 |   .8486625   .0086541   -16.09   0.000     .8318693    .8657948
        vid3 |   1.160412   .0090696    19.04   0.000     1.142771    1.178325
        vid4 |   1.038661   .0065391     6.03   0.000     1.025924    1.051557
        vid5 |   .8948181   .0124975    -7.96   0.000     .8706556    .9196512
        vid6 |   1.080876   .0134682     6.24   0.000     1.054798    1.107598
        vid7 |   .8808432   .0120414    -9.28   0.000     .8575557     .904763
        vid8 |   1.257772   .0122384    23.57   0.000     1.234012    1.281989
 femaleXver2 |   1.919249   .6177648     2.03   0.043     1.021288    3.606739
 femaleXver3 |   1.217535   .3518236     0.68   0.496     .6910587    2.145102
       _cons |   .1917587     .02968   -10.67   0.000     .1415821    .2597181
-------------+----------------------------------------------------------------
4            |
ipv_version2 |   2.331048   .6773259     2.91   0.004     1.318925     4.11986
ipv_version3 |   .4761768   .2140141    -1.65   0.099     .1973341    1.149038
      female |   1.254867   .2608198     1.09   0.275     .8349842    1.885893
        vid2 |   .8486625   .0086541   -16.09   0.000     .8318693    .8657948
        vid3 |   1.160412   .0090696    19.04   0.000     1.142771    1.178325
        vid4 |   1.038661   .0065391     6.03   0.000     1.025924    1.051557
        vid5 |   .8948181   .0124975    -7.96   0.000     .8706556    .9196512
        vid6 |   1.080876   .0134682     6.24   0.000     1.054798    1.107598
        vid7 |   .8808432   .0120414    -9.28   0.000     .8575557     .904763
        vid8 |   1.257772   .0122384    23.57   0.000     1.234012    1.281989
 femaleXver2 |   1.919249   .6177648     2.03   0.043     1.021288    3.606739
 femaleXver3 |   1.217535   .3518236     0.68   0.496     .6910587    2.145102
       _cons |   .0712728   .0160598   -11.72   0.000     .0458272    .1108471
------------------------------------------------------------------------------

. test femaleXver2 femaleXver3

 ( 1)  [0]femaleXver2 = 0
 ( 2)  [1]femaleXver2 = 0
 ( 3)  [2]femaleXver2 = 0
 ( 4)  [3]femaleXver2 = 0
 ( 5)  [4]femaleXver2 = 0
 ( 6)  [0]femaleXver3 = 0
 ( 7)  [1]femaleXver3 = 0
 ( 8)  [2]femaleXver3 = 0
 ( 9)  [3]femaleXver3 = 0
 (10)  [4]femaleXver3 = 0
       Constraint 1 dropped
       Constraint 3 dropped
       Constraint 4 dropped
       Constraint 5 dropped
       Constraint 6 dropped
       Constraint 7 dropped
       Constraint 9 dropped
       Constraint 10 dropped

           chi2(  2) =    4.23
         Prob > chi2 =    0.1209

. 
. ***** sensitivity analysis: try different omnibus variables
. 
. poisson belief_one i.ipv_version i.female i.vid, vce(cluster vid) irr

Iteration 0:   log pseudolikelihood = -1151.0824  
Iteration 1:   log pseudolikelihood = -1151.0823  

Poisson regression                              Number of obs     =      1,334
                                                Wald chi2(6)      =          .
                                                Prob > chi2       =          .
Log pseudolikelihood = -1151.0823               Pseudo R2         =     0.0294

                                    (Std. Err. adjusted for 8 clusters in vid)
------------------------------------------------------------------------------
             |               Robust
  belief_one |        IRR   Std. Err.      z    P>|z|     [95% Conf. Interval]
-------------+----------------------------------------------------------------
 ipv_version |
    2:blame  |   1.579144   .0595492    12.12   0.000     1.466638     1.70028
3:blameless  |   .8264193   .0925478    -1.70   0.089     .6635556    1.029256
             |
      female |
   1:female  |   1.263068   .0414812     7.11   0.000     1.184328    1.347043
             |
         vid |
          2  |   .9315298   .0012931   -51.09   0.000     .9289987    .9340678
          3  |   1.009455   .0006092    15.59   0.000     1.008262    1.010649
          4  |   .9899569   .0011648    -8.58   0.000     .9876766    .9922425
          5  |   .9559877   .0016809   -25.60   0.000     .9526988    .9592879
          6  |   .9014144   .0029643   -31.56   0.000     .8956232    .9072431
          7  |   .9878597   .0015929    -7.58   0.000     .9847426    .9909866
          8  |   .9989074   .0011616    -0.94   0.347     .9966333    1.001187
             |
       _cons |   .4480241   .0217969   -16.50   0.000     .4072766    .4928483
------------------------------------------------------------------------------

. est store belief_one

. 
. poisson norm_one i.ipv_version i.female i.vid, vce(cluster vid) irr

Iteration 0:   log pseudolikelihood = -1206.4256  
Iteration 1:   log pseudolikelihood = -1206.4256  

Poisson regression                              Number of obs     =      1,334
                                                Wald chi2(2)      =          .
                                                Prob > chi2       =          .
Log pseudolikelihood = -1206.4256               Pseudo R2         =     0.0139

                                    (Std. Err. adjusted for 8 clusters in vid)
------------------------------------------------------------------------------
             |               Robust
    norm_one |        IRR   Std. Err.      z    P>|z|     [95% Conf. Interval]
-------------+----------------------------------------------------------------
 ipv_version |
    2:blame  |   1.247635   .0410455     6.73   0.000     1.169726    1.330733
3:blameless  |   .8282397   .0405842    -3.85   0.000     .7523964    .9117283
             |
      female |
   1:female  |   1.205011   .0657905     3.42   0.001     1.082723     1.34111
             |
         vid |
          2  |   .9674707   .0026163   -12.23   0.000     .9623564    .9726123
          3  |   1.090764   .0007563   125.30   0.000     1.089282    1.092247
          4  |   1.013523   .0012292    11.08   0.000     1.011116    1.015935
          5  |     .99037   .0035222    -2.72   0.007     .9834906    .9972976
          6  |   1.090137   .0057608    16.33   0.000     1.078904    1.101487
          7  |   .9615833   .0012246   -30.76   0.000     .9591862    .9639863
          8  |   1.167843   .0010971   165.17   0.000     1.165695    1.169996
             |
       _cons |   .5215301   .0232658   -14.59   0.000     .4778667     .569183
------------------------------------------------------------------------------

. est store norm_one

. 
. poisson belief_all i.ipv_version i.female i.vid, vce(cluster vid) irr

Iteration 0:   log pseudolikelihood = -367.20413  
Iteration 1:   log pseudolikelihood = -367.10121  
Iteration 2:   log pseudolikelihood = -367.10093  
Iteration 3:   log pseudolikelihood = -367.10093  

Poisson regression                              Number of obs     =      1,334
                                                Wald chi2(2)      =          .
                                                Prob > chi2       =          .
Log pseudolikelihood = -367.10093               Pseudo R2         =     0.1279

                                    (Std. Err. adjusted for 8 clusters in vid)
------------------------------------------------------------------------------
             |               Robust
  belief_all |        IRR   Std. Err.      z    P>|z|     [95% Conf. Interval]
-------------+----------------------------------------------------------------
 ipv_version |
    2:blame  |   3.297007   .7390052     5.32   0.000     2.124849     5.11578
3:blameless  |   .2606458   .1053213    -3.33   0.001      .118059    .5754434
             |
      female |
   1:female  |   1.865957   .3139938     3.71   0.000     1.341731    2.595001
             |
         vid |
          2  |   1.451126   .0109478    49.35   0.000     1.429827    1.472743
          3  |   1.730003   .0061486   154.22   0.000     1.717994    1.742096
          4  |   1.241026    .006293    42.58   0.000     1.228753    1.253421
          5  |   1.184847   .0108609    18.50   0.000      1.16375    1.206326
          6  |   1.166243   .0158629    11.31   0.000     1.135563    1.197752
          7  |   1.475334   .0124445    46.10   0.000     1.451144    1.499928
          8  |   1.882036   .0090819   131.04   0.000      1.86432     1.89992
             |
       _cons |   .0296375   .0069503   -15.00   0.000     .0187165    .0469308
------------------------------------------------------------------------------

. est store belief_all

. 
. poisson norm_all i.ipv_version i.female i.vid, vce(cluster vid) irr

Iteration 0:   log pseudolikelihood = -448.44623  
Iteration 1:   log pseudolikelihood = -448.43998  
Iteration 2:   log pseudolikelihood = -448.43998  

Poisson regression                              Number of obs     =      1,334
                                                Wald chi2(2)      =          .
                                                Prob > chi2       =          .
Log pseudolikelihood = -448.43998               Pseudo R2         =     0.0863

                                    (Std. Err. adjusted for 8 clusters in vid)
------------------------------------------------------------------------------
             |               Robust
    norm_all |        IRR   Std. Err.      z    P>|z|     [95% Conf. Interval]
-------------+----------------------------------------------------------------
 ipv_version |
    2:blame  |   2.991626   .6211614     5.28   0.000      1.99145    4.494125
3:blameless  |    .563778    .224113    -1.44   0.149     .2586655     1.22879
             |
      female |
   1:female  |     1.5179   .3520894     1.80   0.072     .9633852    2.391588
             |
         vid |
          2  |   .7151821   .0079336   -30.22   0.000     .6998004    .7309019
          3  |   1.278423   .0048316    64.99   0.000     1.268988    1.287927
          4  |   1.146501   .0040764    38.45   0.000     1.138539    1.154518
          5  |   .8151289   .0075359   -22.11   0.000     .8004918    .8300336
          6  |    1.11886   .0232292     5.41   0.000     1.074246    1.165328
          7  |   1.068367   .0088124     8.02   0.000     1.051234    1.085779
          8  |   1.413813   .0067283    72.77   0.000     1.400687    1.427062
             |
       _cons |   .0552335   .0159251   -10.04   0.000     .0313893    .0971903
------------------------------------------------------------------------------

. est store norm_all

. 
. ***** sensitivity analysis: analyze each item separately 
. 
. foreach x of numlist 1/5 {
  2. display "* analysis for belief `x'"
  3. poisson belief`x'_yes i.ipv_version i.female i.vid, vce(cluster vid) irr
  4. est store belief`x'
  5. }
* analysis for belief 1

Iteration 0:   log pseudolikelihood = -825.50704  
Iteration 1:   log pseudolikelihood = -825.50699  

Poisson regression                              Number of obs     =      1,334
                                                Wald chi2(5)      =          .
                                                Prob > chi2       =          .
Log pseudolikelihood = -825.50699               Pseudo R2         =     0.0313

                                    (Std. Err. adjusted for 8 clusters in vid)
------------------------------------------------------------------------------
             |               Robust
 belief1_yes |        IRR   Std. Err.      z    P>|z|     [95% Conf. Interval]
-------------+----------------------------------------------------------------
 ipv_version |
    2:blame  |   1.813234   .1899546     5.68   0.000     1.476665    2.226516
3:blameless  |   .9316486   .1533884    -0.43   0.667     .6746968    1.286458
             |
      female |
   1:female  |   1.374432   .1175453     3.72   0.000     1.162321    1.625251
             |
         vid |
          2  |   1.061765   .0049978    12.73   0.000     1.052015    1.071606
          3  |   1.198972   .0025815    84.28   0.000     1.193923    1.204042
          4  |   1.261075   .0039351    74.34   0.000     1.253386    1.268811
          5  |   1.092556   .0059657    16.21   0.000     1.080926    1.104312
          6  |   .7856334    .005571   -34.02   0.000     .7747899    .7966286
          7  |   1.203997   .0044193    50.58   0.000     1.195367     1.21269
          8  |   1.284764   .0035292    91.22   0.000     1.277865      1.2917
             |
       _cons |   .1690295   .0155575   -19.31   0.000     .1411296     .202445
------------------------------------------------------------------------------
* analysis for belief 2

Iteration 0:   log pseudolikelihood = -809.48775  
Iteration 1:   log pseudolikelihood = -809.47597  
Iteration 2:   log pseudolikelihood = -809.47597  

Poisson regression                              Number of obs     =      1,334
                                                Wald chi2(2)      =          .
                                                Prob > chi2       =          .
Log pseudolikelihood = -809.47597               Pseudo R2         =     0.1082

                                    (Std. Err. adjusted for 8 clusters in vid)
------------------------------------------------------------------------------
             |               Robust
 belief2_yes |        IRR   Std. Err.      z    P>|z|     [95% Conf. Interval]
-------------+----------------------------------------------------------------
 ipv_version |
    2:blame  |   2.901223   .1243985    24.84   0.000     2.667371    3.155578
3:blameless  |   .5954495    .094154    -3.28   0.001     .4367679    .8117813
             |
      female |
   1:female  |   1.433224   .1262043     4.09   0.000     1.206037    1.703208
             |
         vid |
          2  |   1.021502   .0045932     4.73   0.000     1.012539    1.030544
          3  |   1.187717   .0019939   102.48   0.000     1.183816    1.191631
          4  |   1.145669   .0014325   108.76   0.000     1.142864     1.14848
          5  |   .9319345   .0032758   -20.05   0.000      .925536    .9383772
          6  |   .9537866   .0066823    -6.75   0.000      .940779    .9669741
          7  |   1.178149   .0029799    64.82   0.000     1.172323    1.184004
          8  |   1.189704   .0011934   173.16   0.000     1.187367    1.192045
             |
       _cons |   .1576499   .0099621   -29.23   0.000     .1392853    .1784359
------------------------------------------------------------------------------
* analysis for belief 3

Iteration 0:   log pseudolikelihood = -946.49743  
Iteration 1:   log pseudolikelihood = -946.49326  
Iteration 2:   log pseudolikelihood = -946.49326  

Poisson regression                              Number of obs     =      1,334
                                                Wald chi2(2)      =          .
                                                Prob > chi2       =          .
Log pseudolikelihood = -946.49326               Pseudo R2         =     0.0697

                                    (Std. Err. adjusted for 8 clusters in vid)
------------------------------------------------------------------------------
             |               Robust
 belief3_yes |        IRR   Std. Err.      z    P>|z|     [95% Conf. Interval]
-------------+----------------------------------------------------------------
 ipv_version |
    2:blame  |   2.010802   .1151431    12.20   0.000     1.797329    2.249629
3:blameless  |   .5806531   .0868163    -3.64   0.000     .4331611    .7783663
             |
      female |
   1:female  |   1.287809    .080318     4.06   0.000      1.13963    1.455255
             |
         vid |
          2  |   .9527851   .0025008   -18.43   0.000     .9478962    .9576992
          3  |   1.083977   .0008057   108.49   0.000     1.082399    1.085558
          4  |   1.058736   .0015825    38.19   0.000     1.055639    1.061842
          5  |   .7709183   .0021566   -93.01   0.000     .7667031    .7751567
          6  |   .9354065   .0056462   -11.06   0.000     .9244053    .9465386
          7  |   1.012727   .0026469     4.84   0.000     1.007552    1.017928
          8  |   1.050149   .0018975    27.08   0.000     1.046437    1.053875
             |
       _cons |   .2867142   .0233558   -15.34   0.000      .244405    .3363476
------------------------------------------------------------------------------
* analysis for belief 4

Iteration 0:   log pseudolikelihood = -721.68679  
Iteration 1:   log pseudolikelihood = -721.68654  
Iteration 2:   log pseudolikelihood = -721.68654  

Poisson regression                              Number of obs     =      1,334
                                                Wald chi2(2)      =          .
                                                Prob > chi2       =          .
Log pseudolikelihood = -721.68654               Pseudo R2         =     0.0442

                                    (Std. Err. adjusted for 8 clusters in vid)
------------------------------------------------------------------------------
             |               Robust
 belief4_yes |        IRR   Std. Err.      z    P>|z|     [95% Conf. Interval]
-------------+----------------------------------------------------------------
 ipv_version |
    2:blame  |   2.251908   .3194439     5.72   0.000     1.705311    2.973703
3:blameless  |   .8610761   .1237017    -1.04   0.298     .6497681    1.141103
             |
      female |
   1:female  |   1.291511   .0595988     5.54   0.000     1.179826    1.413768
             |
         vid |
          2  |   1.109633   .0037853    30.50   0.000     1.102238    1.117076
          3  |   1.155123   .0030281    55.01   0.000     1.149204    1.161074
          4  |   1.337302   .0068679    56.60   0.000     1.323908     1.35083
          5  |   1.044049    .007477     6.02   0.000     1.029497    1.058807
          6  |   1.025068   .0061265     4.14   0.000      1.01313    1.037146
          7  |   1.263638   .0065415    45.20   0.000     1.250882    1.276524
          8  |    1.16229   .0038424    45.49   0.000     1.154783    1.169846
             |
       _cons |   .1266261    .011122   -23.53   0.000     .1066005    .1504137
------------------------------------------------------------------------------
* analysis for belief 5

Iteration 0:   log pseudolikelihood = -864.43849  
Iteration 1:   log pseudolikelihood = -864.34766  
Iteration 2:   log pseudolikelihood = -864.34762  

Poisson regression                              Number of obs     =      1,334
                                                Wald chi2(3)      =          .
                                                Prob > chi2       =          .
Log pseudolikelihood = -864.34762               Pseudo R2         =     0.1178

                                    (Std. Err. adjusted for 8 clusters in vid)
------------------------------------------------------------------------------
             |               Robust
 belief5_yes |        IRR   Std. Err.      z    P>|z|     [95% Conf. Interval]
-------------+----------------------------------------------------------------
 ipv_version |
    2:blame  |   2.058389   .1208307    12.30   0.000     1.834681    2.309374
3:blameless  |   .3244069   .0578821    -6.31   0.000      .228673    .4602197
             |
      female |
   1:female  |   1.530815   .0937086     6.96   0.000      1.35774    1.725953
             |
         vid |
          2  |   .8180662   .0021133   -77.74   0.000     .8139347    .8222186
          3  |   .9443546   .0009272   -58.31   0.000      .942539    .9461737
          4  |    .902429   .0020975   -44.17   0.000     .8983273    .9065494
          5  |   .7801529   .0026432   -73.28   0.000     .7749895    .7853506
          6  |   .8499651   .0041035   -33.67   0.000     .8419602     .858046
          7  |   .8085328   .0025501   -67.39   0.000       .80355    .8135464
          8  |    1.05056   .0017809    29.10   0.000     1.047075    1.054056
             |
       _cons |   .2758395   .0164665   -21.57   0.000     .2453823    .3100771
------------------------------------------------------------------------------

. 
. foreach x of numlist 1/5 {
  2. display "* analysis for norm `x'"
  3. poisson norm`x'_yes i.ipv_version i.female i.vid, vce(cluster vid) irr
  4. est store norm`x'
  5. }
* analysis for norm 1

Iteration 0:   log pseudolikelihood = -959.03191  
Iteration 1:   log pseudolikelihood =  -959.0319  

Poisson regression                              Number of obs     =      1,334
                                                Wald chi2(2)      =          .
                                                Prob > chi2       =          .
Log pseudolikelihood =  -959.0319               Pseudo R2         =     0.0151

                                    (Std. Err. adjusted for 8 clusters in vid)
------------------------------------------------------------------------------
             |               Robust
   norm1_yes |        IRR   Std. Err.      z    P>|z|     [95% Conf. Interval]
-------------+----------------------------------------------------------------
 ipv_version |
    2:blame  |    1.36739    .086349     4.96   0.000     1.208204     1.54755
3:blameless  |   .8177122   .0587456    -2.80   0.005     .7103118    .9413519
             |
      female |
   1:female  |   1.175277   .0926266     2.05   0.040     1.007059    1.371594
             |
         vid |
          2  |     .89182   .0037932   -26.92   0.000     .8844164    .8992856
          3  |   1.049647   .0014556    34.94   0.000     1.046798    1.052504
          4  |   .9579081   .0024129   -17.07   0.000     .9531905    .9626489
          5  |   .9170515   .0051344   -15.47   0.000     .9070432    .9271703
          6  |   .8359069   .0061831   -24.23   0.000     .8238757    .8481137
          7  |   1.009164   .0024194     3.81   0.000     1.004433    1.013917
          8  |   1.144925   .0017092    90.66   0.000      1.14158     1.14828
             |
       _cons |   .3177014   .0163137   -22.33   0.000     .2872835    .3513399
------------------------------------------------------------------------------
* analysis for norm 2

Iteration 0:   log pseudolikelihood = -863.92163  
Iteration 1:   log pseudolikelihood = -863.91971  
Iteration 2:   log pseudolikelihood = -863.91971  

Poisson regression                              Number of obs     =      1,334
                                                Wald chi2(2)      =          .
                                                Prob > chi2       =          .
Log pseudolikelihood = -863.91971               Pseudo R2         =     0.0623

                                    (Std. Err. adjusted for 8 clusters in vid)
------------------------------------------------------------------------------
             |               Robust
   norm2_yes |        IRR   Std. Err.      z    P>|z|     [95% Conf. Interval]
-------------+----------------------------------------------------------------
 ipv_version |
    2:blame  |   2.093725   .2033524     7.61   0.000     1.730801     2.53275
3:blameless  |   .6580192   .1155161    -2.38   0.017     .4664542    .9282568
             |
      female |
   1:female  |   1.412607   .1341524     3.64   0.000     1.172693    1.701603
             |
         vid |
          2  |   .9841136   .0050366    -3.13   0.002     .9742914    .9940347
          3  |   1.051757   .0024937    21.28   0.000     1.046881    1.056656
          4  |   1.075407   .0031888    24.52   0.000     1.069175    1.081675
          5  |   .9849312   .0055976    -2.67   0.008     .9740211    .9959636
          6  |   1.104354   .0080242    13.66   0.000     1.088738    1.120193
          7  |   .9060355   .0032274   -27.70   0.000     .8997319    .9123833
          8  |   1.174344   .0029657    63.64   0.000     1.168546    1.180171
             |
       _cons |   .2020447   .0156774   -20.61   0.000       .17354    .2352313
------------------------------------------------------------------------------
* analysis for norm 3

Iteration 0:   log pseudolikelihood =  -1033.364  
Iteration 1:   log pseudolikelihood = -1033.3639  

Poisson regression                              Number of obs     =      1,334
                                                Wald chi2(6)      =          .
                                                Prob > chi2       =          .
Log pseudolikelihood = -1033.3639               Pseudo R2         =     0.0257

                                    (Std. Err. adjusted for 8 clusters in vid)
------------------------------------------------------------------------------
             |               Robust
   norm3_yes |        IRR   Std. Err.      z    P>|z|     [95% Conf. Interval]
-------------+----------------------------------------------------------------
 ipv_version |
    2:blame  |   1.469901   .0646394     8.76   0.000     1.348516    1.602212
3:blameless  |   .7460423   .0494621    -4.42   0.000     .6551327    .8495669
             |
      female |
   1:female  |   1.238193    .119308     2.22   0.027     1.025108    1.495572
             |
         vid |
          2  |   .9529627   .0049119    -9.35   0.000      .943384    .9626386
          3  |   1.078482    .001558    52.30   0.000     1.075433     1.08154
          4  |     1.1039   .0024056    45.36   0.000     1.099195    1.108625
          5  |   1.002867   .0069365     0.41   0.679     .9893636    1.016555
          6  |   1.111183   .0106763    10.97   0.000     1.090453    1.132306
          7  |   .8569006   .0011608  -114.01   0.000     .8546286    .8591787
          8  |   1.108435   .0013341    85.53   0.000     1.105823    1.111053
             |
       _cons |   .3431833   .0203044   -18.08   0.000     .3056081    .3853784
------------------------------------------------------------------------------
* analysis for norm 4

Iteration 0:   log pseudolikelihood = -922.88643  
Iteration 1:   log pseudolikelihood = -922.88643  

Poisson regression                              Number of obs     =      1,334
                                                Wald chi2(2)      =          .
                                                Prob > chi2       =          .
Log pseudolikelihood = -922.88643               Pseudo R2         =     0.0149

                                    (Std. Err. adjusted for 8 clusters in vid)
------------------------------------------------------------------------------
             |               Robust
   norm4_yes |        IRR   Std. Err.      z    P>|z|     [95% Conf. Interval]
-------------+----------------------------------------------------------------
 ipv_version |
    2:blame  |   1.297553   .0987588     3.42   0.001     1.117735    1.506299
3:blameless  |   .8971435   .0296171    -3.29   0.001     .8409331    .9571112
             |
      female |
   1:female  |   1.410365   .1443264     3.36   0.001     1.154053    1.723603
             |
         vid |
          2  |   .9616294   .0049165    -7.65   0.000     .9520413    .9713141
          3  |     1.1886   .0016139   127.24   0.000     1.185441    1.191767
          4  |   1.148112   .0030045    52.78   0.000     1.142239    1.154016
          5  |      .9212   .0057165   -13.23   0.000     .9100637    .9324726
          6  |   1.208485   .0103623    22.08   0.000     1.188344    1.228966
          7  |   1.015877    .002994     5.34   0.000     1.010026    1.021763
          8  |   1.147216   .0019852    79.37   0.000     1.143331    1.151113
             |
       _cons |   .2362185   .0176248   -19.34   0.000     .2040817     .273416
------------------------------------------------------------------------------
* analysis for norm 5

Iteration 0:   log pseudolikelihood =  -879.2029  
Iteration 1:   log pseudolikelihood = -879.19741  
Iteration 2:   log pseudolikelihood = -879.19741  

Poisson regression                              Number of obs     =      1,334
                                                Wald chi2(2)      =          .
                                                Prob > chi2       =          .
Log pseudolikelihood = -879.19741               Pseudo R2         =     0.0701

                                    (Std. Err. adjusted for 8 clusters in vid)
------------------------------------------------------------------------------
             |               Robust
   norm5_yes |        IRR   Std. Err.      z    P>|z|     [95% Conf. Interval]
-------------+----------------------------------------------------------------
 ipv_version |
    2:blame  |   1.942198    .052513    24.55   0.000     1.841954    2.047897
3:blameless  |   .5253366    .058422    -5.79   0.000     .4224513    .6532789
             |
      female |
   1:female  |   1.389253   .1492009     3.06   0.002     1.125551    1.714736
             |
         vid |
          2  |   .8121849   .0044403   -38.05   0.000     .8035285    .8209346
          3  |   1.006312    .001697     3.73   0.000     1.002992    1.009644
          4  |   .9504852   .0017494   -27.59   0.000     .9470626    .9539201
          5  |   .8073145   .0053901   -32.06   0.000     .7968189    .8179482
          6  |   .9928721   .0098926    -0.72   0.473      .973671    1.012452
          7  |   1.007233   .0012856     5.65   0.000     1.004716    1.009756
          8  |   1.018773   .0005328    35.56   0.000      1.01773    1.019818
             |
       _cons |   .2540434   .0148094   -23.51   0.000     .2266143    .2847925
------------------------------------------------------------------------------

. 
. ****** sensitivity analysis: RA/participant sex concordance
. 
. generate femaleXra=female*ra_female

. 
. * beliefs: first show that interaction terms do not vary across logits, then use partia
> l proportional odds
. 
. quietly ologit belief_tot ipv_version2 ipv_version3 female vid2 vid3 vid4 vid5 vid6 vid
> 7 vid8 /*
> */ ra_female femaleXra, vce(cluster vid)

. brant 

Brant Test of Parallel Regression Assumption

    Variable |      chi2   p>chi2    df
-------------+--------------------------
         All |     88.65    0.000    48
-------------+--------------------------
ipv_version2 |     11.31    0.023     4
ipv_version3 |     21.01    0.000     4
      female |      1.95    0.745     4
        vid2 |      3.05    0.549     4
        vid3 |      2.98    0.560     4
        vid4 |      9.62    0.047     4
        vid5 |      4.15    0.386     4
        vid6 |      1.93    0.748     4
        vid7 |      1.27    0.867     4
        vid8 |      5.13    0.274     4
   ra_female |      9.84    0.043     4
   femaleXra |      6.79    0.147     4
----------------------------------------

A significant test statistic provides evidence that the parallel
regression assumption has been violated.

. 
. xi: gologit2 belief_tot ipv_version2 ipv_version3 female vid2 vid3 vid4 vid5 vid6 vid7 
> vid8 /*
> */ ra_female femaleXra, /*
> */ robust cluster(vid) or pl(female vid2 vid3 vid4 vid5 vid6 vid7 vid8 ra_female female
> Xra)

Generalized Ordered Logit Estimates             Number of obs     =      1,334
                                                Wald chi2(2)      =      10.17
                                                Prob > chi2       =     0.0062
Log pseudolikelihood = -1932.0927               Pseudo R2         =     0.0893

 ( 1)  [0]female - [1]female = 0
 ( 2)  [0]vid2 - [1]vid2 = 0
 ( 3)  [0]vid3 - [1]vid3 = 0
 ( 4)  [0]vid4 - [1]vid4 = 0
 ( 5)  [0]vid5 - [1]vid5 = 0
 ( 6)  [0]vid6 - [1]vid6 = 0
 ( 7)  [0]vid7 - [1]vid7 = 0
 ( 8)  [0]vid8 - [1]vid8 = 0
 ( 9)  [0]ra_female - [1]ra_female = 0
 (10)  [0]femaleXra - [1]femaleXra = 0
 (11)  [1]female - [2]female = 0
 (12)  [1]vid2 - [2]vid2 = 0
 (13)  [1]vid3 - [2]vid3 = 0
 (14)  [1]vid4 - [2]vid4 = 0
 (15)  [1]vid5 - [2]vid5 = 0
 (16)  [1]vid6 - [2]vid6 = 0
 (17)  [1]vid7 - [2]vid7 = 0
 (18)  [1]vid8 - [2]vid8 = 0
 (19)  [1]ra_female - [2]ra_female = 0
 (20)  [1]femaleXra - [2]femaleXra = 0
 (21)  [2]female - [3]female = 0
 (22)  [2]vid2 - [3]vid2 = 0
 (23)  [2]vid3 - [3]vid3 = 0
 (24)  [2]vid4 - [3]vid4 = 0
 (25)  [2]vid5 - [3]vid5 = 0
 (26)  [2]vid6 - [3]vid6 = 0
 (27)  [2]vid7 - [3]vid7 = 0
 (28)  [2]vid8 - [3]vid8 = 0
 (29)  [2]ra_female - [3]ra_female = 0
 (30)  [2]femaleXra - [3]femaleXra = 0
 (31)  [3]female - [4]female = 0
 (32)  [3]vid2 - [4]vid2 = 0
 (33)  [3]vid3 - [4]vid3 = 0
 (34)  [3]vid4 - [4]vid4 = 0
 (35)  [3]vid5 - [4]vid5 = 0
 (36)  [3]vid6 - [4]vid6 = 0
 (37)  [3]vid7 - [4]vid7 = 0
 (38)  [3]vid8 - [4]vid8 = 0
 (39)  [3]ra_female - [4]ra_female = 0
 (40)  [3]femaleXra - [4]femaleXra = 0
                                    (Std. Err. adjusted for 8 clusters in vid)
------------------------------------------------------------------------------
             |               Robust
  belief_tot | Odds Ratio   Std. Err.      z    P>|z|     [95% Conf. Interval]
-------------+----------------------------------------------------------------
0            |
ipv_version2 |   3.810873   .4716843    10.81   0.000     2.989981    4.857138
ipv_version3 |   .7019875   .1463243    -1.70   0.090     .4665545    1.056225
      female |   1.767648   .2644356     3.81   0.000     1.318434    2.369918
        vid2 |   .9273229    .012267    -5.70   0.000     .9035891    .9516802
        vid3 |    1.22419   .0186535    13.28   0.000      1.18817    1.261301
        vid4 |   1.072955   .0365567     2.07   0.039     1.003646    1.147052
        vid5 |    .917253   .0067495   -11.74   0.000     .9041192    .9305777
        vid6 |   .8305038   .0184001    -8.38   0.000     .7952121    .8673618
        vid7 |   1.128042   .0129705    10.48   0.000     1.102904    1.153752
        vid8 |   1.175016   .0130155    14.56   0.000     1.149781    1.200804
   ra_female |   .7092595   .1109302    -2.20   0.028     .5220051     .963686
   femaleXra |   1.039359   .1829206     0.22   0.826     .7361387    1.467478
       _cons |   .8845148   .1001096    -1.08   0.278     .7085425    1.104191
-------------+----------------------------------------------------------------
1            |
ipv_version2 |   4.495252   .3947522    17.12   0.000     3.784473    5.339524
ipv_version3 |   .4170289   .0826018    -4.42   0.000      .282857    .6148447
      female |   1.767648   .2644356     3.81   0.000     1.318434    2.369918
        vid2 |   .9273229    .012267    -5.70   0.000     .9035891    .9516802
        vid3 |    1.22419   .0186535    13.28   0.000      1.18817    1.261301
        vid4 |   1.072955   .0365567     2.07   0.039     1.003646    1.147052
        vid5 |    .917253   .0067495   -11.74   0.000     .9041192    .9305777
        vid6 |   .8305038   .0184001    -8.38   0.000     .7952121    .8673618
        vid7 |   1.128042   .0129705    10.48   0.000     1.102904    1.153752
        vid8 |   1.175016   .0130155    14.56   0.000     1.149781    1.200804
   ra_female |   .7092595   .1109302    -2.20   0.028     .5220051     .963686
   femaleXra |   1.039359   .1829206     0.22   0.826     .7361387    1.467478
       _cons |   .4870152   .0422982    -8.28   0.000     .4107844    .5773923
-------------+----------------------------------------------------------------
2            |
ipv_version2 |   5.661253   .7758499    12.65   0.000      4.32772    7.405698
ipv_version3 |   .4549099    .113681    -3.15   0.002     .2787482    .7424014
      female |   1.767648   .2644356     3.81   0.000     1.318434    2.369918
        vid2 |   .9273229    .012267    -5.70   0.000     .9035891    .9516802
        vid3 |    1.22419   .0186535    13.28   0.000      1.18817    1.261301
        vid4 |   1.072955   .0365567     2.07   0.039     1.003646    1.147052
        vid5 |    .917253   .0067495   -11.74   0.000     .9041192    .9305777
        vid6 |   .8305038   .0184001    -8.38   0.000     .7952121    .8673618
        vid7 |   1.128042   .0129705    10.48   0.000     1.102904    1.153752
        vid8 |   1.175016   .0130155    14.56   0.000     1.149781    1.200804
   ra_female |   .7092595   .1109302    -2.20   0.028     .5220051     .963686
   femaleXra |   1.039359   .1829206     0.22   0.826     .7361387    1.467478
       _cons |   .2280407   .0355152    -9.49   0.000     .1680528    .3094418
-------------+----------------------------------------------------------------
3            |
ipv_version2 |   5.118076   .9564184     8.74   0.000     3.548478    7.381954
ipv_version3 |   .3685102   .1163432    -3.16   0.002     .1984791     .684202
      female |   1.767648   .2644356     3.81   0.000     1.318434    2.369918
        vid2 |   .9273229    .012267    -5.70   0.000     .9035891    .9516802
        vid3 |    1.22419   .0186535    13.28   0.000      1.18817    1.261301
        vid4 |   1.072955   .0365567     2.07   0.039     1.003646    1.147052
        vid5 |    .917253   .0067495   -11.74   0.000     .9041192    .9305777
        vid6 |   .8305038   .0184001    -8.38   0.000     .7952121    .8673618
        vid7 |   1.128042   .0129705    10.48   0.000     1.102904    1.153752
        vid8 |   1.175016   .0130155    14.56   0.000     1.149781    1.200804
   ra_female |   .7092595   .1109302    -2.20   0.028     .5220051     .963686
   femaleXra |   1.039359   .1829206     0.22   0.826     .7361387    1.467478
       _cons |   .1228985   .0240696   -10.70   0.000     .0837219    .1804071
-------------+----------------------------------------------------------------
4            |
ipv_version2 |   3.859354   .9689348     5.38   0.000      2.35945    6.312748
ipv_version3 |   .2524276    .102822    -3.38   0.001     .1136095    .5608658
      female |   1.767648   .2644356     3.81   0.000     1.318434    2.369918
        vid2 |   .9273229    .012267    -5.70   0.000     .9035891    .9516802
        vid3 |    1.22419   .0186535    13.28   0.000      1.18817    1.261301
        vid4 |   1.072955   .0365567     2.07   0.039     1.003646    1.147052
        vid5 |    .917253   .0067495   -11.74   0.000     .9041192    .9305777
        vid6 |   .8305038   .0184001    -8.38   0.000     .7952121    .8673618
        vid7 |   1.128042   .0129705    10.48   0.000     1.102904    1.153752
        vid8 |   1.175016   .0130155    14.56   0.000     1.149781    1.200804
   ra_female |   .7092595   .1109302    -2.20   0.028     .5220051     .963686
   femaleXra |   1.039359   .1829206     0.22   0.826     .7361387    1.467478
       _cons |     .05461     .00912   -17.41   0.000     .0393659    .0757573
------------------------------------------------------------------------------

. 
. * norms: first show that interaction terms do not vary across logits, then use partial 
> proportional odds
. 
. quietly ologit norm_tot ipv_version2 ipv_version3 female vid2 vid3 vid4 vid5 vid6 vid7 
> vid8 /*
> */ ra_female femaleXra, vce(cluster vid)

. brant

Brant Test of Parallel Regression Assumption

    Variable |      chi2   p>chi2    df
-------------+--------------------------
         All |     80.40    0.002    48
-------------+--------------------------
ipv_version2 |      9.60    0.048     4
ipv_version3 |      1.97    0.742     4
      female |      5.06    0.281     4
        vid2 |      2.04    0.729     4
        vid3 |      2.45    0.654     4
        vid4 |      0.13    0.998     4
        vid5 |      1.35    0.853     4
        vid6 |      1.27    0.866     4
        vid7 |      0.63    0.959     4
        vid8 |      8.05    0.090     4
   ra_female |      3.42    0.491     4
   femaleXra |      5.87    0.209     4
----------------------------------------

A significant test statistic provides evidence that the parallel
regression assumption has been violated.

. 
. xi: gologit2 norm_tot ipv_version2 ipv_version3 female vid2 vid3 vid4 vid5 vid6 vid7 vi
> d8 /*
> */ ra_female femaleXra, /*
> */ robust cluster(vid) or pl(female vid2 vid3 vid4 vid5 vid6 vid7 vid8 ra_female female
> Xra)

Generalized Ordered Logit Estimates             Number of obs     =      1,334
                                                Wald chi2(2)      =      36.86
                                                Prob > chi2       =     0.0000
Log pseudolikelihood = -2121.6081               Pseudo R2         =     0.0427

 ( 1)  [0]female - [1]female = 0
 ( 2)  [0]vid2 - [1]vid2 = 0
 ( 3)  [0]vid3 - [1]vid3 = 0
 ( 4)  [0]vid4 - [1]vid4 = 0
 ( 5)  [0]vid5 - [1]vid5 = 0
 ( 6)  [0]vid6 - [1]vid6 = 0
 ( 7)  [0]vid7 - [1]vid7 = 0
 ( 8)  [0]vid8 - [1]vid8 = 0
 ( 9)  [0]ra_female - [1]ra_female = 0
 (10)  [0]femaleXra - [1]femaleXra = 0
 (11)  [1]female - [2]female = 0
 (12)  [1]vid2 - [2]vid2 = 0
 (13)  [1]vid3 - [2]vid3 = 0
 (14)  [1]vid4 - [2]vid4 = 0
 (15)  [1]vid5 - [2]vid5 = 0
 (16)  [1]vid6 - [2]vid6 = 0
 (17)  [1]vid7 - [2]vid7 = 0
 (18)  [1]vid8 - [2]vid8 = 0
 (19)  [1]ra_female - [2]ra_female = 0
 (20)  [1]femaleXra - [2]femaleXra = 0
 (21)  [2]female - [3]female = 0
 (22)  [2]vid2 - [3]vid2 = 0
 (23)  [2]vid3 - [3]vid3 = 0
 (24)  [2]vid4 - [3]vid4 = 0
 (25)  [2]vid5 - [3]vid5 = 0
 (26)  [2]vid6 - [3]vid6 = 0
 (27)  [2]vid7 - [3]vid7 = 0
 (28)  [2]vid8 - [3]vid8 = 0
 (29)  [2]ra_female - [3]ra_female = 0
 (30)  [2]femaleXra - [3]femaleXra = 0
 (31)  [3]female - [4]female = 0
 (32)  [3]vid2 - [4]vid2 = 0
 (33)  [3]vid3 - [4]vid3 = 0
 (34)  [3]vid4 - [4]vid4 = 0
 (35)  [3]vid5 - [4]vid5 = 0
 (36)  [3]vid6 - [4]vid6 = 0
 (37)  [3]vid7 - [4]vid7 = 0
 (38)  [3]vid8 - [4]vid8 = 0
 (39)  [3]ra_female - [4]ra_female = 0
 (40)  [3]femaleXra - [4]femaleXra = 0
                                    (Std. Err. adjusted for 8 clusters in vid)
------------------------------------------------------------------------------
             |               Robust
    norm_tot | Odds Ratio   Std. Err.      z    P>|z|     [95% Conf. Interval]
-------------+----------------------------------------------------------------
0            |
ipv_version2 |   2.053053   .2703789     5.46   0.000     1.585989    2.657663
ipv_version3 |    .645422    .073796    -3.83   0.000     .5158452    .8075476
      female |   1.530437   .4268145     1.53   0.127     .8859905    2.643638
        vid2 |   .8529464   .0132893   -10.21   0.000     .8272935    .8793947
        vid3 |   1.140912   .0231332     6.50   0.000     1.096461    1.187165
        vid4 |    1.01826   .0330786     0.56   0.578     .9554481    1.085201
        vid5 |   .8766568   .0133682    -8.63   0.000     .8508432    .9032535
        vid6 |   1.088142   .0297856     3.09   0.002     1.031302    1.148115
        vid7 |   .8724553   .0153606    -7.75   0.000     .8428626    .9030869
        vid8 |   1.250208    .024125    11.57   0.000     1.203807    1.298398
   ra_female |   .8968082   .1935299    -0.50   0.614     .5875049     1.36895
   femaleXra |   1.136981   .3419267     0.43   0.669     .6306225    2.049919
       _cons |   1.223186   .2060442     1.20   0.232     .8792425    1.701674
-------------+----------------------------------------------------------------
1            |
ipv_version2 |   2.211936   .3200235     5.49   0.000     1.665791     2.93714
ipv_version3 |   .6277537   .0623513    -4.69   0.000     .5167067    .7626662
      female |   1.530437   .4268145     1.53   0.127     .8859905    2.643638
        vid2 |   .8529464   .0132893   -10.21   0.000     .8272935    .8793947
        vid3 |   1.140912   .0231332     6.50   0.000     1.096461    1.187165
        vid4 |    1.01826   .0330786     0.56   0.578     .9554481    1.085201
        vid5 |   .8766568   .0133682    -8.63   0.000     .8508432    .9032535
        vid6 |   1.088142   .0297856     3.09   0.002     1.031302    1.148115
        vid7 |   .8724553   .0153606    -7.75   0.000     .8428626    .9030869
        vid8 |   1.250208    .024125    11.57   0.000     1.203807    1.298398
   ra_female |   .8968082   .1935299    -0.50   0.614     .5875049     1.36895
   femaleXra |   1.136981   .3419267     0.43   0.669     .6306225    2.049919
       _cons |   .6462647   .0904231    -3.12   0.002     .4912622    .8501734
-------------+----------------------------------------------------------------
2            |
ipv_version2 |   2.842168   .3624425     8.19   0.000     2.213614    3.649201
ipv_version3 |   .5623937   .0730212    -4.43   0.000     .4360341    .7253715
      female |   1.530437   .4268145     1.53   0.127     .8859905    2.643638
        vid2 |   .8529464   .0132893   -10.21   0.000     .8272935    .8793947
        vid3 |   1.140912   .0231332     6.50   0.000     1.096461    1.187165
        vid4 |    1.01826   .0330786     0.56   0.578     .9554481    1.085201
        vid5 |   .8766568   .0133682    -8.63   0.000     .8508432    .9032535
        vid6 |   1.088142   .0297856     3.09   0.002     1.031302    1.148115
        vid7 |   .8724553   .0153606    -7.75   0.000     .8428626    .9030869
        vid8 |   1.250208    .024125    11.57   0.000     1.203807    1.298398
   ra_female |   .8968082   .1935299    -0.50   0.614     .5875049     1.36895
   femaleXra |   1.136981   .3419267     0.43   0.669     .6306225    2.049919
       _cons |   .3370928   .0479619    -7.64   0.000     .2550584     .445512
-------------+----------------------------------------------------------------
3            |
ipv_version2 |   3.054801   .3050054    11.18   0.000     2.511857    3.715102
ipv_version3 |   .4874997   .0813769    -4.30   0.000      .351468    .6761809
      female |   1.530437   .4268145     1.53   0.127     .8859905    2.643638
        vid2 |   .8529464   .0132893   -10.21   0.000     .8272935    .8793947
        vid3 |   1.140912   .0231332     6.50   0.000     1.096461    1.187165
        vid4 |    1.01826   .0330786     0.56   0.578     .9554481    1.085201
        vid5 |   .8766568   .0133682    -8.63   0.000     .8508432    .9032535
        vid6 |   1.088142   .0297856     3.09   0.002     1.031302    1.148115
        vid7 |   .8724553   .0153606    -7.75   0.000     .8428626    .9030869
        vid8 |   1.250208    .024125    11.57   0.000     1.203807    1.298398
   ra_female |   .8968082   .1935299    -0.50   0.614     .5875049     1.36895
   femaleXra |   1.136981   .3419267     0.43   0.669     .6306225    2.049919
       _cons |   .1752759   .0342241    -8.92   0.000     .1195413     .256996
-------------+----------------------------------------------------------------
4            |
ipv_version2 |   3.507682   .8110166     5.43   0.000     2.229526    5.518586
ipv_version3 |   .5346785   .2251258    -1.49   0.137     .2342585    1.220366
      female |   1.530437   .4268145     1.53   0.127     .8859905    2.643638
        vid2 |   .8529464   .0132893   -10.21   0.000     .8272935    .8793947
        vid3 |   1.140912   .0231332     6.50   0.000     1.096461    1.187165
        vid4 |    1.01826   .0330786     0.56   0.578     .9554481    1.085201
        vid5 |   .8766568   .0133682    -8.63   0.000     .8508432    .9032535
        vid6 |   1.088142   .0297856     3.09   0.002     1.031302    1.148115
        vid7 |   .8724553   .0153606    -7.75   0.000     .8428626    .9030869
        vid8 |   1.250208    .024125    11.57   0.000     1.203807    1.298398
   ra_female |   .8968082   .1935299    -0.50   0.614     .5875049     1.36895
   femaleXra |   1.136981   .3419267     0.43   0.669     .6306225    2.049919
       _cons |   .0649029   .0130656   -13.59   0.000      .043743    .0962985
------------------------------------------------------------------------------

. 
. ***** sensitivity analysis: use PCA to generate dependent variables (response to review
> er #4)
. 
. quietly pca belief1_yes belief2_yes belief3_yes belief4_yes belief5_yes, correlation

. /* quietly predict belief_index */
. summarize belief_index, detail

          Belief index (first principal component)
-------------------------------------------------------------
      Percentiles      Smallest
 1%    -1.536719      -1.536719
 5%    -1.536719      -1.536719
10%    -1.536719      -1.536719       Obs               1,334
25%    -1.536719      -1.536719       Sum of Wgt.       1,334

50%    -.5861626                      Mean           1.10e-08
                        Largest       Std. Dev.      1.703602
75%     1.428955       3.294702
90%     2.425261       3.294702       Variance       2.902261
95%     3.294702       3.294702       Skewness       .6937545
99%     3.294702       3.294702       Kurtosis       2.036643

. regress belief_index i.ipv_version i.female i.vid, vce(cluster vid) 

Linear regression                               Number of obs     =      1,334
                                                F(2, 7)           =          .
                                                Prob > F          =          .
                                                R-squared         =     0.2525
                                                Root MSE          =     1.4784

                                    (Std. Err. adjusted for 8 clusters in vid)
------------------------------------------------------------------------------
             |               Robust
belief_index |      Coef.   Std. Err.      t    P>|t|     [95% Conf. Interval]
-------------+----------------------------------------------------------------
 ipv_version |
    2:blame  |   1.441591   .0737104    19.56   0.000     1.267294    1.615889
3:blameless  |  -.4694494   .1108521    -4.23   0.004    -.7315728   -.2073259
             |
      female |
   1:female  |   .4860248   .0485124    10.02   0.000     .3713113    .6007384
             |
         vid |
          2  |  -.0492876   .0030771   -16.02   0.000    -.0565638   -.0420114
          3  |   .1470363   .0014833    99.13   0.000     .1435288    .1505438
          4  |   .1514321   .0025753    58.80   0.000     .1453425    .1575218
          5  |  -.1657246   .0046799   -35.41   0.000    -.1767909   -.1546583
          6  |  -.1393968   .0054075   -25.78   0.000    -.1521836     -.12661
          7  |   .0751916   .0022578    33.30   0.000     .0698526    .0805305
          8  |   .1891315   .0019935    94.88   0.000     .1844177    .1938453
             |
       _cons |  -.6289157   .0416614   -15.10   0.000    -.7274292   -.5304022
------------------------------------------------------------------------------

. 
. quietly pca norm1_yes norm2_yes norm3_yes norm4_yes norm5_yes, correlation

. /* quietly predict norm_index */
. summarize norm_index, detail

           Norm index (first principal component)
-------------------------------------------------------------
      Percentiles      Smallest
 1%     -1.67604       -1.67604
 5%     -1.67604       -1.67604
10%     -1.67604       -1.67604       Obs               1,334
25%     -1.67604       -1.67604       Sum of Wgt.       1,334

50%    -.6826862                      Mean           4.30e-09
                        Largest       Std. Dev.      1.686321
75%     1.256221       2.988991
90%     2.988991       2.988991       Variance        2.84368
95%     2.988991       2.988991       Skewness       .5181802
99%     2.988991       2.988991       Kurtosis       1.822784

. regress norm_index i.ipv_version i.female i.vid, vce(cluster vid) 

Linear regression                               Number of obs     =      1,334
                                                F(2, 7)           =          .
                                                Prob > F          =          .
                                                R-squared         =     0.1320
                                                Root MSE          =      1.577

                                    (Std. Err. adjusted for 8 clusters in vid)
------------------------------------------------------------------------------
             |               Robust
  norm_index |      Coef.   Std. Err.      t    P>|t|     [95% Conf. Interval]
-------------+----------------------------------------------------------------
 ipv_version |
    2:blame  |   .9325033   .0825746    11.29   0.000     .7372454    1.127761
3:blameless  |  -.4061258   .0675009    -6.02   0.001    -.5657402   -.2465114
             |
      female |
   1:female  |   .4474894   .1105245     4.05   0.005     .1861404    .7088385
             |
         vid |
          2  |  -.1384962   .0063864   -21.69   0.000    -.1535976   -.1233948
          3  |   .1174068   .0017429    67.36   0.000     .1132856     .121528
          4  |   .0704708   .0032885    21.43   0.000     .0626947     .078247
          5  |  -.1213505   .0084537   -14.35   0.000    -.1413403   -.1013608
          6  |   .0770676   .0108231     7.12   0.000     .0514751    .1026601
          7  |   -.075839   .0024397   -31.09   0.000    -.0816079   -.0700701
          8  |   .1868316   .0019235    97.13   0.000     .1822832      .19138
             |
       _cons |  -.4530394   .0554119    -8.18   0.000    -.5840678   -.3220111
------------------------------------------------------------------------------

. 
. log close
      name:  <unnamed>
       log:  /Users/drdrtsai/Documents/Stata/Data/HopeNet/IPV_RCT/ipvrct_replication_2017
> 0408.log
  log type:  text
 closed on:   9 Apr 2017, 13:32:21
-----------------------------------------------------------------------------------------
